Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,475,822 T→C M129T (ATG→ACG)  yjgL → SopA‑central‑domain‑like hexapeptide repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,475,8220TC99.0% 354.5 / ‑4.2 95M129T (ATG→ACG) yjgLSopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (50/44);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

AATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGTTACCAAAGTGAAG  >  NC_000913/4475674‑4475968
                                                                                                                                                    |                                                                                                                                                  
aaTATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACgg                                                                                                                                                  <  2:183935/151‑1 (MQ=255)
   atgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACgggg                                                                                                                                                >  2:760079/1‑150 (MQ=255)
   atgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACgggg                                                                                                                                                >  2:37911/1‑150 (MQ=255)
    tgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGAc                                                                                                                                              >  1:1174499/1‑151 (MQ=255)
        ggCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGt                                                                                                                                          >  1:911451/1‑151 (MQ=255)
             atatTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAgg                                                                                                                                     <  1:442969/151‑1 (MQ=255)
             atatTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAgg                                                                                                                                     <  2:744214/151‑1 (MQ=255)
               atTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGa                                                                                                                                    <  1:981447/150‑1 (MQ=255)
               atTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGAc                                                                                                                                   >  2:1636724/1‑151 (MQ=255)
                  aaGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACAt                                                                                                                                 <  2:588876/150‑1 (MQ=255)
                      ccTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAg                                                                                                                             >  2:142678/1‑150 (MQ=255)
                       cTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGcc                                                                                                                           >  1:249480/1‑151 (MQ=255)
                          aTCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCAtt                                                                                                                        >  1:1124350/1‑151 (MQ=255)
                           tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCAtt                                                                                                                        >  2:1430819/1‑150 (MQ=255)
                              tATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCAtt                                                                                                                        <  2:1460757/147‑1 (MQ=255)
                                 tCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATAttt                                                                                                                  <  1:693025/150‑1 (MQ=255)
                                    aaaGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTcc                                                                                                                <  1:1086493/149‑1 (MQ=255)
                                    aaaGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTc                                                                                                                 <  2:127644/148‑1 (MQ=255)
                                      aGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCcaca                                                                                                             >  2:256620/1‑150 (MQ=255)
                                      aGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCcaca                                                                                                             >  1:957734/1‑150 (MQ=255)
                                      aGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACaa                                                                                                            <  2:240270/151‑1 (MQ=255)
                                          aCTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTg                                                                                                        >  2:985767/1‑151 (MQ=255)
                                            tAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGAt                                                                                                      >  1:1007894/1‑151 (MQ=255)
                                             aaTGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATg                                                                                                     >  1:876147/1‑151 (MQ=255)
                                             aaTGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATg                                                                                                     <  1:1012869/151‑1 (MQ=255)
                                               tGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGt                                                                                                   >  2:1252150/1‑151 (MQ=255)
                                                ggCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGtt                                                                                                  <  2:605027/151‑1 (MQ=255)
                                                ggCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGtt                                                                                                  <  2:346632/151‑1 (MQ=255)
                                                 gCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTc                                                                                                 <  1:1092554/151‑1 (MQ=255)
                                                  cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCg                                                                                                >  1:1500161/1‑151 (MQ=255)
                                                  cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCg                                                                                                <  1:1311465/151‑1 (MQ=255)
                                                      gAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACt                                                                                             <  2:126713/150‑1 (MQ=255)
                                                      gAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTa                                                                                            <  1:83098/151‑1 (MQ=255)
                                                          aaaTGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCg                                                                                        >  2:1212820/1‑151 (MQ=255)
                                                            aTGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCCCTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGtt                                                                                      <  1:215491/151‑1 (MQ=255)
                                                            aTGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGtt                                                                                      >  2:1446406/1‑151 (MQ=255)
                                                             tGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGtt                                                                                      >  2:1247420/1‑150 (MQ=255)
                                                                   tttACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAg                                                                                <  1:699441/150‑1 (MQ=255)
                                                                     tACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAggg                                                                              >  1:1307053/1‑150 (MQ=255)
                                                                      aCATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGaa                                                                            <  2:1126441/151‑1 (MQ=255)
                                                                       cATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGa                                                                             >  2:1500245/1‑149 (MQ=255)
                                                                       cATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGa                                                                             >  2:808390/1‑149 (MQ=255)
                                                                       cATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTACCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAg                                                                           <  2:1513653/151‑1 (MQ=255)
                                                                        aTGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGt                                                                          >  1:179712/1‑151 (MQ=255)
                                                                        aTGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGt                                                                          >  2:108435/1‑151 (MQ=255)
                                                                         tGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTc                                                                         >  1:1364323/1‑151 (MQ=255)
                                                                         tGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTc                                                                         >  2:785651/1‑151 (MQ=255)
                                                                          gAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTcc                                                                        <  1:297747/151‑1 (MQ=255)
                                                                          gAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTcc                                                                        >  2:662560/1‑151 (MQ=255)
                                                                              aatCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTg                                                                    >  2:1008843/1‑151 (MQ=255)
                                                                              aatCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTg                                                                    <  1:1688806/151‑1 (MQ=255)
                                                                               atCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTga                                                                   <  2:957807/151‑1 (MQ=255)
                                                                               atCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTg                                                                    >  2:1678650/1‑150 (MQ=255)
                                                                                tCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTga                                                                   >  1:969552/1‑150 (MQ=255)
                                                                                 ccTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGCCCTTTgaaa                                                                 <  2:464189/151‑1 (MQ=255)
                                                                                  cTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTgaa                                                                  >  2:1546212/1‑149 (MQ=255)
                                                                                   tGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGtt                                                                                                  <  2:1467463/116‑1 (MQ=255)
                                                                                   tGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGtt                                                                                                  >  1:1467365/1‑116 (MQ=255)
                                                                                   tGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTgaa                                                                  >  1:1382075/1‑148 (MQ=255)
                                                                                    gTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTgaaagaa                                                              <  2:1218421/151‑1 (MQ=255)
                                                                                    gTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTgaaaga                                                               >  1:695541/1‑150 (MQ=255)
                                                                                      aGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTgaaagaaa                                                             <  1:168980/150‑1 (MQ=255)
                                                                                      aGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGaaagaaaa                                                            <  1:1636619/151‑1 (MQ=255)
                                                                                        gggAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTgaaagaaa                                                             >  1:540604/1‑148 (MQ=255)
                                                                                           aaCGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATg                                                        >  2:193622/1‑150 (MQ=255)
                                                                                            aCGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATgag                                                      <  1:1279404/151‑1 (MQ=255)
                                                                                              gAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATgag                                                      >  1:348617/1‑149 (MQ=255)
                                                                                              gAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGagag                                                    <  1:214485/151‑1 (MQ=255)
                                                                                                 cGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTc                                                 <  1:735384/151‑1 (MQ=255)
                                                                                                   gATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCAt                                               <  1:1705128/151‑1 (MQ=255)
                                                                                                      tCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTa                                             <  2:1268848/150‑1 (MQ=255)
                                                                                                          ccAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCt                                        >  2:887501/1‑151 (MQ=255)
                                                                                                           cAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTg                                       >  1:946926/1‑151 (MQ=255)
                                                                                                           cAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTg                                       <  1:1173462/151‑1 (MQ=255)
                                                                                                            aaGGTGATTATAGAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTg                                       <  1:243855/150‑1 (MQ=255)
                                                                                                                     ataGAAATTTCACTTTCCACTATCCCTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAAc                              <  1:1617691/150‑1 (MQ=255)
                                                                                                                        gAAATTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTTCACCTGATTACTAACAAgg                          >  1:1418227/1‑151 (MQ=255)
                                                                                                                         aaaTTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGt                         >  2:1663713/1‑151 (MQ=255)
                                                                                                                         aaaTTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGt                         <  2:1424946/151‑1 (MQ=255)
                                                                                                                          aaTTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTTCACCTGATTACTAACAAGGtt                        <  2:1418320/151‑1 (MQ=255)
                                                                                                                           aTTTCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTc                       >  1:801767/1‑151 (MQ=255)
                                                                                                                             ttCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCt                      <  1:1338433/150‑1 (MQ=255)
                                                                                                                             ttCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCt                      <  2:92062/150‑1 (MQ=255)
                                                                                                                             ttCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTg                     >  2:93360/1‑151 (MQ=255)
                                                                                                                              tCACTTTCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTgg                    >  2:1241260/1‑151 (MQ=255)
                                                                                                                                   ttCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCgtg                <  2:911519/150‑1 (MQ=255)
                                                                                                                                    tCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGtt              >  2:617464/1‑151 (MQ=255)
                                                                                                                                    tCCACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGtt              >  1:891486/1‑151 (MQ=255)
                                                                                                                                     ccACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGtt              >  2:655928/1‑150 (MQ=255)
                                                                                                                                     ccACTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGtt              >  2:769065/1‑150 (MQ=255)
                                                                                                                                       aCTATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGTTAcc           >  2:1195923/1‑151 (MQ=255)
                                                                                                                                         tATCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGTTACCa          >  1:924915/1‑150 (MQ=255)
                                                                                                                                          aTCACTACTACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGTTACCaaa        <  2:577340/151‑1 (MQ=255)
                                                                                                                                             actactACGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGTTACCAAAGTg     <  2:249503/151‑1 (MQ=255)
                                                                                                                                                actaCGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGTTACCAAAGTGAAg  <  2:109791/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGAAGTCCTTTGAAAGAAAAGATGAGAGTTCATTACACCTGATTACTAACAAGGTTCTGGCGTGTTACCAAAGTGAAG  >  NC_000913/4475674‑4475968

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