Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,839,567 G→C 5.4% L143V (CTT→GTT)  nlpA ← lipoprotein‑28

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,839,5670GC 5.4% 312.4 / 7.1 92L143V (CTT→GTT) nlpAlipoprotein‑28
Reads supporting (aligned to +/- strand):  ref base G (47/40);  new base C (3/2);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

AGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGTT  >  NC_000913/3839419‑3839713
                                                                                                                                                    |                                                                                                                                                  
aGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAgg                                                                                                                                                   <  2:228736/150‑1 (MQ=255)
 gTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGt                                                                                                                                                  >  2:158323/1‑150 (MQ=255)
   aCGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTg                                                                                                                                                <  1:1021354/150‑1 (MQ=255)
      cgcgGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTggtc                                                                                                                                             >  1:1637558/1‑150 (MQ=255)
          ggAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAACCTTggtcgcgc                                                                                                                                         <  1:321918/150‑4 (MQ=255)
           gAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTggtcgggtc                                                                                                                                        >  2:980686/1‑150 (MQ=255)
              cTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCCTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGtt                                                                                                                                     >  2:1256445/1‑150 (MQ=255)
               tGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTg                                                                                                                                    >  1:819194/1‑150 (MQ=255)
               tGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTg                                                                                                                                    >  1:423010/1‑150 (MQ=255)
               tGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTg                                                                                                                                    <  1:387618/150‑1 (MQ=255)
               tGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTg                                                                                                                                    >  1:181405/1‑150 (MQ=255)
                gcgcCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTgg                                                                                                                                   <  1:176435/150‑1 (MQ=255)
                gcgcCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTgg                                                                                                                                   <  1:159861/150‑1 (MQ=255)
                 cgcCCCCTCGAGTTCCATAATCTGCAAATGAGGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTggg                                                                                                                                  >  2:989339/1‑150 (MQ=255)
                       cTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGc                                                                                                                            >  1:444203/1‑150 (MQ=255)
                        tCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGcc                                                                                                                           <  2:1460900/150‑1 (MQ=255)
                         cGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCa                                                                                                                          >  2:237672/1‑150 (MQ=255)
                         cGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCa                                                                                                                          >  2:475493/1‑150 (MQ=255)
                          gAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCAc                                                                                                                         <  2:106808/150‑1 (MQ=255)
                            gTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTg                                                                                                                       >  1:804208/1‑150 (MQ=255)
                               ccATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcg                                                                                                                    >  1:1617276/1‑150 (MQ=255)
                               ccATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcg                                                                                                                    >  1:238325/1‑150 (MQ=255)
                                cATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgc                                                                                                                   <  2:181405/150‑1 (MQ=255)
                                cATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgc                                                                                                                   <  1:1429628/150‑1 (MQ=255)
                                 aTAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGGAACAGTAAAAGTGCGCGACCAAGGTTggtcggg                                                                                                                                          >  1:908876/1‑126 (MQ=255)
                                 aTAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTggtcggg                                                                                                                                          <  2:908876/126‑1 (MQ=255)
                                 aTAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgcg                                                                                                                  <  1:158323/150‑1 (MQ=255)
                                  tAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCgcgc                                                                                                                 >  2:958003/1‑150 (MQ=255)
                                         cAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCttt                                                                                                          >  2:579541/1‑150 (MQ=255)
                                         cAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCttt                                                                                                          >  1:1476921/1‑150 (MQ=255)
                                             tGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTAtt                                                                                                      <  2:819194/150‑1 (MQ=255)
                                             tGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTAGCGCCTTCTTTTAtt                                                                                                      <  1:1630779/150‑1 (MQ=255)
                                               aCGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTg                                                                                                    <  1:813859/150‑1 (MQ=255)
                                                 gAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTgcg                                                                                                  <  1:1571309/150‑1 (MQ=255)
                                                          cGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGttt                                                                                         <  1:231050/150‑1 (MQ=255)
                                                           ggTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCct                                                                                                               <  1:220434/127‑1 (MQ=255)
                                                           ggTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCct                                                                                                               >  2:220434/1‑127 (MQ=255)
                                                           ggTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGGCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGtttt                                                                                        <  1:1576078/150‑1 (MQ=255)
                                                               atatCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAtt                                                                                    >  2:1292423/1‑150 (MQ=255)
                                                                tatCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAttt                                                                                   >  1:1088257/1‑150 (MQ=255)
                                                                tatCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAttt                                                                                   >  1:1082535/1‑150 (MQ=255)
                                                                   ccAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTGGAtttttt                                                                                >  2:760565/1‑150 (MQ=255)
                                                                     aGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTg                                                                              <  1:67268/150‑1 (MQ=255)
                                                                       cgcgGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGa                                                                            <  1:1014085/150‑1 (MQ=255)
                                                                              ggtaATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAAccc                                                                     <  2:1363947/150‑1 (MQ=255)
                                                                                taataaCCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGc                                                                   >  1:1250981/1‑150 (MQ=255)
                                                                                 aataaCCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTTTTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGcc                                                                  >  2:1305484/1‑150 (MQ=255)
                                                                                       ccTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAggg                                                            <  1:1592454/150‑1 (MQ=255)
                                                                                       ccTTTTCCCTCATTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCTCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAggg                                                            <  2:756067/150‑1 (MQ=255)
                                                                                        cTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGa                                                           >  1:423266/1‑150 (MQ=255)
                                                                                           ttCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGa                                                        >  1:702958/1‑150 (MQ=255)
                                                                                            tCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGAc                                                       >  1:706012/1‑150 (MQ=255)
                                                                                             cccTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCAt                                                                >  2:259550/1‑140 (MQ=255)
                                                                                             cccTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCAt                                                                <  1:259550/140‑1 (MQ=255)
                                                                                             cccTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACa                                                      >  2:1344236/1‑150 (MQ=255)
                                                                                               ctctTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACaaa                                                    >  1:494822/1‑150 (MQ=255)
                                                                                                 ctTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACg                                                  <  1:1248858/150‑1 (MQ=255)
                                                                                                 ctTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACg                                                  <  1:763051/150‑1 (MQ=255)
                                                                                                     cAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTAtt                                              <  2:201931/150‑1 (MQ=255)
                                                                                                       gagTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTc                                            >  1:1340046/1‑150 (MQ=255)
                                                                                                         gTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTccc                                          >  1:224620/1‑150 (MQ=255)
                                                                                                           aaTTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGGCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCAc                                        <  2:106792/150‑1 (MQ=255)
                                                                                                            aTTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCAcc                                       <  2:238325/150‑1 (MQ=255)
                                                                                                             ttAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGGCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCg                                      <  2:496629/150‑1 (MQ=255)
                                                                                                               aGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCt                                    >  2:50008/1‑150 (MQ=255)
                                                                                                                gCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTa                                   <  2:1402146/150‑1 (MQ=255)
                                                                                                                 cTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                                  <  2:804208/150‑1 (MQ=255)
                                                                                                                      gtttttGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTAGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                             <  1:970363/149‑1 (MQ=255)
                                                                                                                      ctTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                             >  1:595075/1‑150 (MQ=255)
                                                                                                                      ctTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                             >  2:234066/1‑150 (MQ=255)
                                                                                                                        ttttGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTta                           <  1:1344236/150‑1 (MQ=255)
                                                                                                                            gCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAc                       >  1:764295/1‑150 (MQ=255)
                                                                                                                             cAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAcc                      >  2:256730/1‑150 (MQ=255)
                                                                                                                             cAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAcc                      <  2:1617276/150‑1 (MQ=255)
                                                                                                                             cAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAcc                      <  1:475493/150‑1 (MQ=255)
                                                                                                                             cAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAcc                      >  1:1509082/1‑150 (MQ=255)
                                                                                                                              aaCAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCa                     <  1:1256445/150‑1 (MQ=255)
                                                                                                                                cAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATg                   <  1:237672/150‑1 (MQ=255)
                                                                                                                                   tAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                             <  1:794879/137‑1 (MQ=255)
                                                                                                                                   tAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                             >  2:794879/1‑137 (MQ=255)
                                                                                                                                   tAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATgcgc                >  1:281703/1‑150 (MQ=255)
                                                                                                                                       acgacccgACCAACCTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAcc                                 <  2:387214/123‑1 (MQ=255)
                                                                                                                                       acgacccgACCAACCTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAcc                                 >  1:387214/7‑129 (MQ=255)
                                                                                                                                       aGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGa            >  2:1205651/1‑150 (MQ=255)
                                                                                                                                         tGCGCGCCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATACCCCGCCATAGGGAAGACAAACGGATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGAtt          >  1:259695/1‑150 (MQ=255)
                                                                                                                                          gcgcGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTa         <  1:234066/150‑1 (MQ=255)
                                                                                                                                          gcgcGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTa         >  1:752366/1‑150 (MQ=255)
                                                                                                                                          acccgACCAACCTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTa         >  1:513744/4‑150 (MQ=255)
                                                                                                                                          acccgACCAACCTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTa         >  1:1631024/4‑150 (MQ=255)
                                                                                                                                               aCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTg    >  2:469677/1‑150 (MQ=255)
                                                                                                                                                ccAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGt   <  2:1188566/150‑1 (MQ=255)
                                                                                                                                                 cAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGtt  >  1:1259938/1‑150 (MQ=255)
                                                                                                                                                 cAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGtt  >  1:793055/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGTT  >  NC_000913/3839419‑3839713

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