Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 405,463 G→A 5.3% A61V (GCG→GTG)  proC ← pyrroline‑5‑carboxylate reductase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913405,4630GA 5.3% 248.6 / 2.5 75A61V (GCG→GTG) proCpyrroline‑5‑carboxylate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (32/38);  new base A (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCA  >  NC_000913/405320‑405607
                                                                                                                                               |                                                                                                                                                
gcgcgcGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCAc                                                                                                                                              >  6:381336/1‑148 (MQ=255)
 cgcgcgGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCAct                                                                                                                                             >  7:209920/1‑148 (MQ=255)
  gcgcgGGCAATCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCActt                                                                                                                                            <  1:1165/148‑1 (MQ=255)
             cTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGa                                                                                                                                  <  1:45098/147‑1 (MQ=255)
              tGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCTCCACTTCTTTCGCCTAt                                                                                                                                 >  2:45098/1‑147 (MQ=255)
              tGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGAtt                                                                                                                                >  5:270226/1‑148 (MQ=255)
                gTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCt                                                                                                                              <  5:140864/148‑1 (MQ=255)
                 tCGAGCGTGACACCTGCAGCAATAGAAACGCCCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGGTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTGATTGCGCCGAttaag                                                                                                                             >  6:335981/1‑145 (MQ=255)
                     gCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTgcggc                                                                                                                         <  4:246340/148‑1 (MQ=255)
                       gTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGt                                                                                                                       <  3:143526/148‑1 (MQ=255)
                        tGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGtt                                                                                                                      <  2:40310/148‑1 (MQ=255)
                         gACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  5:381336/148‑1 (MQ=255)
                         gACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  8:209920/148‑1 (MQ=255)
                          acacCTGCAGCAATAGAAACGACCAGATAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGAAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  5:335981/147‑1 (MQ=255)
                          acacCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  6:270226/147‑1 (MQ=255)
                          acacCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  3:271202/147‑1 (MQ=255)
                                gcagcaATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     >  2:174498/1‑141 (MQ=255)
                                gcagcaATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  1:174498/141‑1 (MQ=255)
                                  agcaATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAAc                                                                                                            <  1:187010/148‑1 (MQ=255)
                                  agcaATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAAc                                                                                                            <  1:181541/148‑1 (MQ=255)
                                     aaTAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTgg                                                                                                         >  8:391867/1‑148 (MQ=255)
                                        aGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCACCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCa                                                                                                      >  2:240718/1‑148 (MQ=255)
                                           aaCGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  2:251872/130‑1 (MQ=255)
                                           aaCGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATg                                                                                                    >  7:285021/1‑147 (MQ=255)
                                           aaCGACCAGAGAGTCTTTAGTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     >  1:251872/1‑130 (MQ=255)
                                                       gTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAAc                                                                                                            <  6:362892/127‑1 (MQ=255)
                                                       gTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAAc                                                                                                            >  5:362892/1‑127 (MQ=255)
                                                        tCTTTATTCATGCTGGAGGTGATTTCTCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  4:272694/117‑1 (MQ=255)
                                                         cTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACtt                                                                                      >  1:208317/1‑147 (MQ=255)
                                                         cTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACtt                                                                                      >  6:213729/1‑147 (MQ=255)
                                                         cTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCACCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACtt                                                                                      >  5:379529/1‑147 (MQ=255)
                                                         cTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGACAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGa                                                                                                                    >  3:272694/1‑117 (MQ=255)
                                                            tATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCACCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTAt                                                                                   <  6:379529/147‑1 (MQ=255)
                                                             aTTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATc                                                                                  <  8:297588/147‑1 (MQ=255)
                                                              ttCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGa                                                                                                              >  6:67928/1‑118 (MQ=255)
                                                              ttCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGa                                                                                                              <  5:67928/118‑1 (MQ=255)
                                                                cAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGAt                                                                                                                   >  3:231223/1‑111 (MQ=255)
                                                                cAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGAt                                                                                                                   <  4:231223/111‑1 (MQ=255)
                                                                 aGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCACCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCgggg                                                                             <  1:240718/148‑1 (MQ=255)
                                                                   gCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAAc                                                                                                            >  1:144205/1‑115 (MQ=255)
                                                                   gCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAAc                                                                                                            <  2:144205/115‑1 (MQ=255)
                                                                    cTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAgg                                                                          <  1:50757/148‑1 (MQ=255)
                                                                       gAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTgcg                                                                                                                                                  <  8:301682/73‑1 (MQ=255)
                                                                       gAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTgcg                                                                                                                                                  >  7:301682/1‑73 (MQ=255)
                                                                          gTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGgtgt                                                                    <  8:285021/148‑1 (MQ=255)
                                                                           tGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     <  3:329994/98‑1 (MQ=255)
                                                                           tGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTg                                                                                                                     >  4:329994/1‑98 (MQ=255)
                                                                              tttCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATAc                                                                >  6:209545/1‑148 (MQ=255)
                                                                                            ttAATCGTGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGc                                                                                                   >  5:212904/1‑99 (MQ=255)
                                                                                            ttAATCGTGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGc                                                                                                   <  6:212904/99‑1 (MQ=255)
                                                                                                tCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAg                                              >  6:383662/1‑148 (MQ=255)
                                                                                                          aGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCg                                    >  2:49938/1‑148 (MQ=255)
                                                                                                          aGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCg                                    <  7:391867/148‑1 (MQ=255)
                                                                                                          aGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCg                                    >  8:351185/1‑148 (MQ=255)
                                                                                                             tttAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGCTGCCGAACTGGTCATGCAGGGCGGCGACTGTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGCCCGCTg                                 >  3:270509/1‑148 (MQ=255)
                                                                                                             tttAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGgtgt                                                                    >  7:298508/1‑113 (MQ=255)
                                                                                                             tttAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGgtgt                                                                    <  8:298508/113‑1 (MQ=255)
                                                                                                             tttAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTg                                                       <  4:112095/126‑1 (MQ=255)
                                                                                                             tttAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTg                                                       >  3:112095/1‑126 (MQ=255)
                                                                                                             tttAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGAGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTg                                 >  7:33464/1‑148 (MQ=255)
                                                                                                               tAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTgg                                >  3:85373/1‑147 (MQ=255)
                                                                                                                 aCGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCa                              <  5:383662/147‑1 (MQ=255)
                                                                                                                   ggCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATc                           <  5:213729/148‑1 (MQ=255)
                                                                                                                     cagcaAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAg                         <  1:49938/148‑1 (MQ=255)
                                                                                                                     cagcaAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAg                         >  5:71101/1‑148 (MQ=255)
                                                                                                                     cagcaAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAg                         <  7:351185/148‑1 (MQ=255)
                                                                                                                      agcaAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGa                        <  6:71101/148‑1 (MQ=255)
                                                                                                                        caAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGa                                                <  7:314274/122‑1 (MQ=255)
                                                                                                                        caAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGa                                                >  8:314274/1‑122 (MQ=255)
                                                                                                                           aaaTGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGAccgcc                   >  7:175907/1‑148 (MQ=255)
                                                                                                                            aaTGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGAccgcc                   >  7:285296/1‑147 (MQ=255)
                                                                                                                              tgatgTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCgag                <  5:209545/148‑1 (MQ=255)
                                                                                                                                 tgTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAAt             <  8:285296/148‑1 (MQ=255)
                                                                                                                                  gTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATg            >  1:196181/1‑148 (MQ=255)
                                                                                                                                   tCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGtata                                                                 <  6:20992/94‑1 (MQ=255)
                                                                                                                                   tCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGtata                                                                 >  5:20992/1‑94 (MQ=255)
                                                                                                                                      gCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCtt        <  2:208317/148‑1 (MQ=255)
                                                                                                                                             gcgcCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCa  <  8:175907/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCA  >  NC_000913/405320‑405607

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: