Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,424 A→T 9.1% intergenic (‑442/+161) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4240AT 9.1% 165.7 / 6.5 66intergenic (‑442/+161)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base T (3/3);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTG  >  NC_000913/987277‑987573
                                                                                                                                                   |                                                                                                                                                      
cTTAAATTTAAATTAACTAAACTGAAATTTTGATTCTTATTGAAACTTTAAAGACCTGGCACAAAATCTTCTTTCAGACATCCAGACTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTCATCCGCAATAAGAATg                                                                                                                                                      <  2:199385/150‑1 (MQ=255)
  tataTTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCg                                                                                                                                                    >  2:941724/1‑150 (MQ=255)
  tataTTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCg                                                                                                                                                    <  1:144829/150‑1 (MQ=255)
      tttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                                <  2:364732/150‑1 (MQ=255)
       ttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTa                                                                                                                                               <  1:1053096/150‑1 (MQ=255)
          aaTTAACTAAACGGAAATTTTGTTTCTGATGGCAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGtt                                                                                                                                            >  2:168227/1‑150 (MQ=255)
          aaTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCAATTAGCTAAGAAAAAACTAATCAGCAATAAGAATGCGGATTAGtt                                                                                                                                            >  2:1143474/1‑150 (MQ=255)
           aTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGttt                                                                                                                                           <  1:1126922/150‑1 (MQ=255)
             tAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGttttt                                                                                                                                         <  1:1144026/150‑1 (MQ=255)
             tAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGttttt                                                                                                                                         <  1:573025/150‑1 (MQ=255)
              aaCTAAACGGAAATTTTGTTTCTGATTGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGGGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTg                                                                                                                                        <  1:56228/150‑1 (MQ=255)
              aaCTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTg                                                                                                                                        <  2:1138008/150‑1 (MQ=255)
                cTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGttttttct                                                                                                                                      <  1:356958/150‑4 (MQ=39)
                 tAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTTTGCTATTAGCTAAGAAAAAACTAATCCTCAATAAGAATGAGGATTAGTTTTTgcgc                                                                                                                                     >  1:1111614/1‑150 (MQ=255)
                  acaCGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGcgcg                                                                                                                                    <  1:836214/148‑1 (MQ=255)
                  aaaCGGAAATTTTGTTTCTGATGGCAACTTTATCGACCTGGCACAACATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCACAATAAGAATGCGGATTAGTTTTTGcgcg                                                                                                                                    <  1:630008/150‑1 (MQ=255)
                          aTTTTGTTTCTGATGGAAAATTTATAGACCTGGCACAAAATCTTCTTTAAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGtt                                                                                                                            <  2:1146976/150‑1 (MQ=255)
                           ttttGTTTCTGATGGAAACTTTATCGACCTTGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATATTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATACGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTg                                                                                                                           <  1:1144916/150‑1 (MQ=255)
                               gTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAAt                                                                                                                       >  1:1045631/1‑150 (MQ=255)
                                    tGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGACAAAAAATAGTATGAAAATTGTGCTATTAGCTAATAAAAAACTAATCCGCAATAAGAATGCGGATTTGTTTTTTCGCGTAAATGTTGCAATGTAat                                                                                                                  >  2:668501/1‑150 (MQ=255)
                                     gATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAata                                                                                                                 >  1:1121923/1‑150 (MQ=255)
                                       tGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATAtt                                                                                                               >  2:313980/1‑150 (MQ=255)
                                        ggAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTc                                                                                                              >  2:173442/1‑150 (MQ=255)
                                          aaaCTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCg                                                                                                            <  2:92477/150‑1 (MQ=255)
                                                tATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAg                                                                                                      <  1:747644/150‑1 (MQ=255)
                                                 aTCGACCTGGCACAAAATCTTCTTTCAGACAGCCAGAATGCCAAAAAATAGTATGAAAATTGTGGTATTAGCTAAGAAAAAAATAATCCGCAATAAGATGGCGGATTTGTTTTTGGGCGTAAATGTTGCAATGTATTTTTCCGTAACAgg                                                                                                     >  2:528091/1‑150 (MQ=255)
                                                    gACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGgatg                                                                                                  >  2:1297968/1‑150 (MQ=255)
                                                     aCCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGgatga                                                                                                 >  2:824074/1‑150 (MQ=255)
                                                      ccTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGgatgat                                                                                                >  1:496948/1‑150 (MQ=255)
                                                      ccTGGCACAAAATCATCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAATAAAAAACTAATCCGCAATTAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGgatgat                                                                                                <  2:319258/150‑1 (MQ=255)
                                                            acaAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTAt                                                                                          <  1:17501/150‑1 (MQ=255)
                                                              aaaaTCTTCTTTCATACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATgg                                                                                        <  2:1257882/150‑1 (MQ=255)
                                                               aaaTCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGc                                                                                       <  1:74934/150‑1 (MQ=255)
                                                                    ttcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGt                                                                                  <  1:454445/150‑1 (MQ=255)
                                                                    ttcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGt                                                                                  <  1:440641/150‑1 (MQ=255)
                                                                      cttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTg                                                                                <  1:779526/150‑1 (MQ=255)
                                                                        ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCa                                                                              >  1:3669/1‑150 (MQ=255)
                                                                           aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAAt                                                                           <  2:660603/150‑1 (MQ=255)
                                                                               aTCCAGAATGCCAAAAAAAAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGAGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGAt                                                                       <  1:26322/150‑1 (MQ=255)
                                                                                    gAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCtt                                                                  <  1:1000937/150‑1 (MQ=255)
                                                                                         ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCTGAAATGTTGCAATTTATTATTCCGTAAAAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGATTGACtt                                                             >  2:106413/1‑150 (MQ=255)
                                                                                           aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCAGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaa                                                           >  2:183071/1‑150 (MQ=255)
                                                                                             aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaaca                                                         <  2:947748/150‑1 (MQ=255)
                                                                                              aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaa                                                        <  2:335977/150‑1 (MQ=255)
                                                                                                 tAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaaca                                                     <  1:731924/150‑1 (MQ=255)
                                                                                                   gTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGa                                                   <  1:554533/150‑1 (MQ=255)
                                                                                                    tATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGAt                                                  >  2:948949/1‑150 (MQ=255)
                                                                                                     aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAAAAGATTGATTGACTTAACAAACAGATa                                                 >  2:1274546/1‑150 (MQ=255)
                                                                                                      tGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAGATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATaa                                                <  1:777228/150‑1 (MQ=255)
                                                                                                       gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAAc                                               >  2:998036/1‑150 (MQ=255)
                                                                                                             tgtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACa                                         >  2:706839/1‑150 (MQ=255)
                                                                                                              gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGAGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAg                                        <  1:631275/150‑1 (MQ=255)
                                                                                                                  tATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaaa                                    <  2:124225/150‑1 (MQ=255)
                                                                                                                   aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAg                                   <  2:566644/150‑1 (MQ=255)
                                                                                                                       gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCAGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaa                              <  1:630041/150‑1 (MQ=255)
                                                                                                                             aaaaaaCTAATCCTCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAag                        <  1:332276/150‑1 (MQ=255)
                                                                                                                               aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAATATAAAAAAagag                      >  1:406301/1‑150 (MQ=255)
                                                                                                                                 aaCTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAGTATTCCGTAACAGGATGATCGTTATGGCGGGGTTGAAAATAGATTGATTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAATAGcc                    >  1:140010/1‑150 (MQ=255)
                                                                                                                                 aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc                    >  1:603945/1‑150 (MQ=255)
                                                                                                                                   cTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAAGAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAg                  >  1:191445/1‑150 (MQ=255)
                                                                                                                                      aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg               >  2:1081803/1‑150 (MQ=255)
                                                                                                                                      aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTAGCAATGTAATATTCCGTAACAGGTTGCTCGTTAGGGCGGGGTTGCAAATAGATTGCGAGACTTAACAAACAGATAGCTGGACAGAAAAGTTAAATAAAGAGCCCGCgg               >  2:858114/1‑150 (MQ=255)
                                                                                                                                        ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAAGAGCCAGCggg              >  2:799861/1‑150 (MQ=255)
                                                                                                                                        ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc             >  1:878359/1‑150 (MQ=255)
                                                                                                                                             aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg        >  1:740862/1‑150 (MQ=255)
                                                                                                                                             aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg        >  1:1124859/1‑150 (MQ=255)
                                                                                                                                                  aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtt   <  2:849378/150‑1 (MQ=255)
                                                                                                                                                   aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCTTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCttc  >  2:583014/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTG  >  NC_000913/987277‑987573

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