Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,105,388 G→T 5.6% intergenic (+95/‑26) csgC → / → ymdA inhibitor of CsgA amyloid formation/uncharacterized protein YmdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,105,3880GT 5.6% 295.9 / 5.4 89intergenic (+95/‑26)csgC/ymdAinhibitor of CsgA amyloid formation/uncharacterized protein YmdA
Reads supporting (aligned to +/- strand):  ref base G (48/36);  new base T (2/3);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

TGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGG  >  NC_000913/1105239‑1105534
                                                                                                                                                     |                                                                                                                                                  
tGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATg                                                                                                                                                    >  2:277318/1‑150 (MQ=255)
 ggACAGTCACTTCATTTATCACAACAATGGCCGCCGTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGCGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATgc                                                                                                                                                   >  2:1955095/1‑149 (MQ=255)
          cTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAATTCTTAATTTGTTGAAATATCGAGCATAAGATTAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCAtt                                                                                                                                          <  1:823325/150‑1 (MQ=255)
           ttCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTg                                                                                                                                         <  1:1140111/150‑1 (MQ=255)
             cATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGACGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGATTAGTCATCCTTGAg                                                                                                                                       <  2:35637/150‑1 (MQ=255)
              aTTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAgg                                                                                                                                      >  1:1814747/1‑150 (MQ=255)
              aTTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTAGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGATTAGTCATCCTTGAgg                                                                                                                                      >  1:1839210/1‑150 (MQ=255)
              aTTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTAGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGATTAGTCATCCTTGAgg                                                                                                                                      <  2:1839210/150‑1 (MQ=255)
               tttATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAATGATTAGTCATCCTTGAggg                                                                                                                                     <  1:1623595/150‑2 (MQ=255)
               tttATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGa                                                                                                                                     >  1:835504/1‑150 (MQ=255)
               tttATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGa                                                                                                                                     <  2:421217/150‑1 (MQ=255)
                  aTCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAAt                                                                                                                                  >  1:1018971/1‑150 (MQ=255)
                  aTCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAAt                                                                                                                                  >  2:1397144/1‑150 (MQ=255)
                   tCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAATTATGGATAACACAACCCTCAAGTATGACTAATCATTGAGGAAATa                                                                                                                                 <  2:259101/150‑1 (MQ=255)
                   tCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATa                                                                                                                                 >  2:844782/1‑150 (MQ=255)
                   tCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATa                                                                                                                                 >  1:887148/1‑150 (MQ=255)
                    cacaACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTTTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTTAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGTAAATAg                                                                                                                                <  1:230975/150‑1 (MQ=255)
                    cacaACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAg                                                                                                                                >  2:1954413/1‑150 (MQ=255)
                    cacaACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAg                                                                                                                                <  2:771782/150‑1 (MQ=255)
                    cacaACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAg                                                                                                                                >  1:862943/1‑150 (MQ=255)
                     acaacaATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGa                                                                                                                               >  1:1173249/1‑150 (MQ=255)
                      caacaaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGaa                                                                                                                              >  2:446068/1‑150 (MQ=255)
                      caacaaTGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGaa                                                                                                                              >  2:938980/1‑150 (MQ=255)
                             ggCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGt                                                                                                                       >  2:749703/1‑150 (MQ=255)
                              gccgccCTCTTCAGAAAAGGCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCGGAAAGGATGGATCCCACAACCCTCAAGGATGACTAATCATTGAGGAAAGAGAATAAATggt                                                                                                                      >  1:344677/1‑148 (MQ=255)
                                      ctTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCtt                                                                                                              >  2:1369920/1‑150 (MQ=255)
                                          aGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCtttttt                                                                                                          >  2:1164337/1‑150 (MQ=255)
                                           gAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTa                                                                                                         <  2:998652/150‑1 (MQ=255)
                                           gAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTa                                                                                                         >  1:1176934/1‑150 (MQ=255)
                                             aaaGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTaaa                                                                                                       >  1:1100835/1‑150 (MQ=255)
                                               aGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAAct                                                                                                     <  1:844782/150‑1 (MQ=255)
                                                   ttAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAActctct                                                                                                 >  2:1127840/1‑150 (MQ=255)
                                                   ttAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGCATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAActctct                                                                                                 <  2:308634/150‑1 (MQ=255)
                                                     aaTTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTa                                                                                               >  2:962478/1‑150 (MQ=255)
                                                     aaTTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTa                                                                                               <  2:1018971/150‑1 (MQ=255)
                                                      aTTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTAt                                                                                              <  2:997331/150‑1 (MQ=255)
                                                      aTTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTAt                                                                                              >  1:1239458/1‑150 (MQ=255)
                                                         tgttgAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCt                                                                                           >  2:219999/1‑150 (MQ=255)
                                                           ttgAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCg                                                                                         <  2:631538/150‑1 (MQ=255)
                                                             gAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGc                                                                                       <  2:1173249/150‑1 (MQ=255)
                                                             gAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGc                                                                                       >  2:1389868/1‑150 (MQ=255)
                                                             gAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGc                                                                                       >  1:887921/1‑150 (MQ=255)
                                                             gAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTAGGc                                                                                       >  1:881230/1‑150 (MQ=255)
                                                               aaTATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAg                                                                                     <  2:1636797/150‑1 (MQ=255)
                                                                   tCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAgtttg                                                                                 >  2:483994/1‑150 (MQ=255)
                                                                       gCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGtttt                                                                             >  2:196676/1‑150 (MQ=255)
                                                                          tAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTc                                                                          >  1:811442/1‑150 (MQ=255)
                                                                           aaGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTcc                                                                         <  1:196676/150‑1 (MQ=255)
                                                                                   tCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTg                                                                 <  2:771937/150‑1 (MQ=255)
                                                                                   tCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTg                                                                 <  1:1727120/150‑1 (MQ=255)
                                                                                   tCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTg                                                                 <  1:1954413/150‑1 (MQ=255)
                                                                                    cTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGc                                                                <  1:277318/150‑1 (MQ=255)
                                                                                    cTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGc                                                                >  1:1788112/1‑150 (MQ=255)
                                                                                        agagaATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCAtt                                                            >  2:425631/1‑150 (MQ=255)
                                                                                          agaATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGc                                                          <  2:376656/150‑1 (MQ=255)
                                                                                          agaATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGc                                                          <  2:369437/150‑1 (MQ=255)
                                                                                             aTGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGTAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTgg                                                       <  2:1676492/150‑1 (MQ=255)
                                                                                             aTGGTCTGCTGCGAATCAGCCAACCTGAAAGTAGGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAAAAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTgg                                                       <  2:666487/150‑1 (MQ=255)
                                                                                               ggTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGaa                                                     >  2:692527/1‑150 (MQ=255)
                                                                                               ggTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACATTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGGTGGac                                                     >  2:803451/1‑149 (MQ=255)
                                                                                                gTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAg                                                    >  1:1208676/1‑150 (MQ=255)
                                                                                                  ctgctgCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTGTTTTACTTTTATTGCCATTGCTGGAAGcc                                                  >  2:330073/1‑149 (MQ=255)
                                                                                                  ctgctgCGAATCAGCCAACCTGAAAGTATGGAAAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTATTAAACTCTCTTATGCTCGTCAGATTGTTTATTCCTTTTATTCCCATTGCTGGAAGCa                                                  >  2:1130476/1‑150 (MQ=255)
                                                                                                    gctgcGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCAcc                                                <  2:1100835/150‑1 (MQ=255)
                                                                                                     ctgcGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCg                                               >  2:145720/1‑150 (MQ=255)
                                                                                                         gAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCa                                           <  1:1127840/150‑1 (MQ=255)
                                                                                                                ccAACCTGAAAGTATGGATAACCCAACCCTCAAGGATGCGTAATCATTGAGGAAATAGAATAAATGTTCAGCCCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCGAGGGGGc                                    >  1:311285/1‑150 (MQ=255)
                                                                                                                 cAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCg                                   >  1:1228747/1‑150 (MQ=255)
                                                                                                                 cAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCg                                   >  1:1798879/1‑150 (MQ=255)
                                                                                                                  aaCCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGt                                  <  1:1397144/150‑1 (MQ=255)
                                                                                                                    ccTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGa                                >  1:680888/1‑150 (MQ=255)
                                                                                                                    ccTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGa                                >  1:1731608/1‑150 (MQ=255)
                                                                                                                     cTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGAt                               <  2:835504/150‑1 (MQ=255)
                                                                                                                      tGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATc                              >  1:1490207/1‑150 (MQ=255)
                                                                                                                        aaaGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCa                            >  1:1977294/1‑150 (MQ=255)
                                                                                                                          aGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCAtt                          >  2:305242/1‑150 (MQ=255)
                                                                                                                              tGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTa                      <  2:687748/150‑1 (MQ=255)
                                                                                                                              tGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTa                      <  1:938980/150‑1 (MQ=255)
                                                                                                                                gATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATg                    <  1:219999/150‑1 (MQ=255)
                                                                                                                                gATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATg                    <  1:1998957/150‑1 (MQ=255)
                                                                                                                                 aTAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATgg                   <  2:1208676/150‑1 (MQ=255)
                                                                                                                                   aaCACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGcc                 <  2:881230/150‑1 (MQ=255)
                                                                                                                                   aaCACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGcc                 <  2:887921/150‑1 (MQ=255)
                                                                                                                                         aCCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAgggg                                      >  2:1871819/1‑123 (MQ=255)
                                                                                                                                         aCCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAgggg                                      <  1:1871819/123‑1 (MQ=255)
                                                                                                                                           ccTCAATGATTAGTCATCCTTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATtgtg         >  1:695848/1‑150 (MQ=255)
                                                                                                                                              cAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGaa      >  1:294289/1‑150 (MQ=255)
                                                                                                                                              cAAGGATGACTAATAATTGAGGAAATAGAATAAATGTTCAGACCTTATTTAAACTCTCTTATGCTAGGAAATTTGATTTTTCCTTTTATTCCCATTGATGGAATCACAGTTAAAGGGGGAGTGATACATTTTTATGGCCAAATTGTGGaa      >  2:1754100/1‑150 (MQ=255)
                                                                                                                                                 ggATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCg   <  2:1977294/150‑1 (MQ=255)
                                                                                                                                                 ggATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCg   >  2:931500/1‑150 (MQ=255)
                                                                                                                                                  gATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCgg  >  1:1486205/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGG  >  NC_000913/1105239‑1105534

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: