Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,302,452 A→C 8.3% L101R (CTG→CGG)  napA ← nitrate reductase, periplasmic, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,302,4520AC 8.3% 305.9 / 20.5 96L101R (CTG→CGG) napAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base A (44/44);  new base C (5/3);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

CAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCA  >  NC_000913/2302309‑2302597
                                                                                                                                               |                                                                                                                                                 
cAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGc                                                                                                                                              >  3:677242/1‑149 (MQ=255)
cAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGc                                                                                                                                              >  3:223006/1‑149 (MQ=255)
   aaCCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGc                                                                                                                                           <  1:168042/149‑1 (MQ=255)
       gAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGtat                                                                                                                                       <  4:181883/149‑3 (MQ=255)
           aTACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCCGCTGCTGCGtatgaaa                                                                                                                                   <  2:54615/149‑7 (MQ=255)
           aTACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACg                                                                                                                                   >  2:135305/1‑149 (MQ=255)
            tACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACgg                                                                                                                                  >  3:667931/1‑149 (MQ=255)
              ccGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTc                                                                                                                                <  4:779967/149‑1 (MQ=255)
               cgatcgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCt                                                                                                                               >  3:532641/1‑149 (MQ=255)
               cgatcgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCt                                                                                                                               >  1:612657/1‑149 (MQ=255)
                gatcgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCtt                                                                                                                              >  3:851327/1‑149 (MQ=255)
                gatcgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCtt                                                                                                                              >  3:610129/1‑149 (MQ=255)
                gatcgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCtt                                                                                                                              >  2:884867/1‑149 (MQ=255)
                gatcgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCtt                                                                                                                              <  2:85776/149‑1 (MQ=255)
                 atcgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCttt                                                                                                                             <  4:97235/149‑1 (MQ=255)
                 atcgatTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAAAGGTCttt                                                                                                                             >  2:399537/1‑149 (MQ=255)
                     atTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGACCTTCTCTTCCATCTCATCCAAGGCCAGATTCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACACAGCAGCGGCTGCGTCAAACGGTCTTTACCg                                                                                                                         <  2:63795/149‑1 (MQ=255)
                      ttCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGt                                                                                                                        >  3:450138/1‑149 (MQ=255)
                       tCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTCTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTa                                                                                                                       <  2:372456/149‑1 (MQ=255)
                        ccGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTAc                                                                                                                      <  1:135305/149‑1 (MQ=255)
                             ccTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCAGATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGAt                                                                                                                 >  4:509475/1‑149 (MQ=255)
                                ttttCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGAtttt                                                                                                              <  4:532641/149‑1 (MQ=255)
                                 tttctttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTg                                                                                                             <  1:619180/149‑1 (MQ=255)
                                  ttctttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                            <  3:140175/149‑1 (MQ=255)
                                  ttctttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                            <  3:518997/149‑1 (MQ=255)
                                  ttctttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                            <  4:450138/149‑1 (MQ=255)
                                     tttcAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCa                                                                                                         >  3:695015/1‑149 (MQ=255)
                                         aGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGaa                                                                                                     >  3:21939/1‑149 (MQ=255)
                                          gggCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTa                                                                                                                            >  2:453581/1‑125 (MQ=255)
                                          gggCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTa                                                                                                                            <  1:453581/125‑1 (MQ=255)
                                          gggCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGaaa                                                                                                    <  2:821523/149‑1 (MQ=255)
                                          gggCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGaaa                                                                                                    <  4:54040/149‑1 (MQ=255)
                                            gCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATa                                                                                                  <  3:273354/149‑1 (MQ=255)
                                              ggTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGc                                                                                                <  2:587110/149‑1 (MQ=255)
                                               gTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGcc                                                                                               <  4:201969/149‑1 (MQ=255)
                                               gTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGcc                                                                                               <  2:46038/149‑1 (MQ=255)
                                                ttttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGccc                                                                                              >  1:242078/1‑149 (MQ=255)
                                                 tttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACGGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCt                                                                                             >  3:600074/1‑149 (MQ=255)
                                                 tttGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCt                                                                                             >  3:659370/1‑149 (MQ=255)
                                                    gAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTaa                                                                                          <  2:597221/149‑1 (MQ=255)
                                                    gAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTaa                                                                                          <  1:140059/149‑1 (MQ=255)
                                                     gacttCTCTTCCATCACCCCGAAGGCCTTATCCCCAGTGAATGGGGTAAATTCGCCTTTTTTGGCATTTTTTCCGGTTTTCCTACGCCGCAGCGGCCGCCGCAAAACGGCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAAt                                                                                         <  3:74575/146‑1 (MQ=255)
                                                     aaTTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAAt                                                                                         <  3:9991/149‑1 (MQ=255)
                                                       tttCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGc                                                                                       >  2:888634/1‑149 (MQ=255)
                                                          ctctTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGt                                                                                    >  2:622472/1‑149 (MQ=255)
                                                             ttCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCa                                                                                 >  2:680818/1‑149 (MQ=255)
                                                             ttCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCa                                                                                 <  3:514145/149‑1 (MQ=255)
                                                              tCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAACCGGTCTTTCCCGTACATGATTTTGGGCAGGAAATAGCCCTTCATGCCGTTCAg                                                                                <  4:837184/149‑1 (MQ=255)
                                                              tCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAg                                                                                <  2:270978/149‑1 (MQ=255)
                                                              tCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAg                                                                                <  3:244386/149‑1 (MQ=255)
                                                                cATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGc                                                                              >  3:509075/1‑149 (MQ=255)
                                                                          aaGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAAc                                                                    >  1:666878/1‑149 (MQ=255)
                                                                              ccTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGt                                                                <  4:495519/149‑1 (MQ=255)
                                                                               cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTg                                                               >  3:330316/1‑149 (MQ=255)
                                                                               cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTg                                                               >  2:262588/1‑149 (MQ=255)
                                                                               cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTg                                                               >  4:228820/1‑149 (MQ=255)
                                                                               cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTg                                                               >  2:449930/1‑149 (MQ=255)
                                                                               cTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTg                                                               >  3:654077/1‑149 (MQ=255)
                                                                                tGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGc                                                              <  4:100781/149‑1 (MQ=255)
                                                                                tGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGc                                                              <  1:233862/149‑1 (MQ=255)
                                                                                 gATCCCATGTGATTGGGGTAAGTTCGGCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCg                                                             <  4:21939/149‑1 (MQ=255)
                                                                                           gATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGc                                                   <  3:259276/149‑1 (MQ=255)
                                                                                            aTTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGcc                                                  <  3:474865/149‑1 (MQ=255)
                                                                                             ttGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGccc                                                 <  4:667931/149‑1 (MQ=255)
                                                                                                   tAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAg                                           <  4:362120/149‑1 (MQ=255)
                                                                                                       ttCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGcca                                       >  1:664393/1‑149 (MQ=255)
                                                                                                       ttCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGcca                                       <  2:460852/149‑1 (MQ=255)
                                                                                                       ttCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGcca                                       >  1:292875/1‑149 (MQ=255)
                                                                                                           ccTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCcaccac                                   >  4:493645/1‑149 (MQ=255)
                                                                                                            cttcttTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACcaca                                  >  1:251248/1‑149 (MQ=255)
                                                                                                            cttcttTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACcaca                                  >  2:290025/1‑149 (MQ=255)
                                                                                                             ttcttTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACa                                            >  1:326174/1‑138 (MQ=255)
                                                                                                             ttcttTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACa                                            <  2:326174/138‑1 (MQ=255)
                                                                                                                tttGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTc                              >  2:493905/1‑149 (MQ=255)
                                                                                                                tttGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTc                              >  2:638329/1‑149 (MQ=255)
                                                                                                                 ttGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTcc                             <  1:449930/149‑1 (MQ=255)
                                                                                                                    tCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCctg                          >  3:124133/1‑149 (MQ=255)
                                                                                                                           tACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGtt                   >  4:780935/1‑149 (MQ=255)
                                                                                                                            aCCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACg                                                                          >  4:206972/1‑93 (MQ=255)
                                                                                                                            aCCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACg                                                                          <  3:206972/93‑1 (MQ=255)
                                                                                                                             ccGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTcc                 >  3:766939/1‑149 (MQ=255)
                                                                                                                             ccGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTcc                 >  4:751206/1‑149 (MQ=255)
                                                                                                                                tttttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGAc              <  4:223006/149‑1 (MQ=255)
                                                                                                                                tttttCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGAc              >  3:623122/1‑149 (MQ=255)
                                                                                                                                   gtttgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAg           >  4:266911/6‑149 (MQ=255)
                                                                                                                                    tttgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGa          >  1:807983/5‑149 (MQ=255)
                                                                                                                                    tttgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGa          >  1:612126/5‑149 (MQ=255)
                                                                                                                                     cATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGaa         <  4:851327/149‑1 (MQ=255)
                                                                                                                                     cATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGaa         <  3:751206/149‑1 (MQ=255)
                                                                                                                                      tgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCctgctg                       <  1:156072/132‑1 (MQ=255)
                                                                                                                                      tgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCctgctg                       >  2:156072/3‑134 (MQ=255)
                                                                                                                                      aTACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAg                                           <  3:338846/114‑1 (MQ=255)
                                                                                                                                      aTACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAg                                           >  4:338846/1‑114 (MQ=255)
                                                                                                                                       tACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCg                                                                  >  1:581694/1‑90 (MQ=255)
                                                                                                                                       tACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCg                                                                  <  2:581694/90‑1 (MQ=255)
                                                                                                                                        aCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGAGAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgc      >  3:580348/1‑149 (MQ=255)
                                                                                                                                          gcagcaGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgccg    >  4:184507/1‑149 (MQ=255)
                                                                                                                                          gcagcaGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGACACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAccccg    >  2:212306/1‑149 (MQ=255)
                                                                                                                                             gcagcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCGGCGTTCTGGTCGGAACgcagca  <  1:809763/148‑1 (MQ=255)
                                                                                                                                             gcagcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCGGCGTTCTGGTCGGAACgcagca  >  2:809763/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCA  >  NC_000913/2302309‑2302597

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: