Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,420 A→C 6.8% intergenic (‑438/+165) 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,4200AC 6.8% 182.5 / 3.9 68intergenic (‑438/+165)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base C (2/3);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGG  >  NC_000913/987271‑987567
                                                                                                                                                     |                                                                                                                                                    
taTTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGAAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATa                                                                                                                                                      <  2:368842/149‑1 (MQ=255)
 aTTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATaa                                                                                                                                                     >  2:1100886/1‑150 (MQ=255)
  tttCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAg                                                                                                                                                    <  2:222543/150‑1 (MQ=255)
   ttCCTTATATTTTAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCAACCTGGCACAAAATCTTCTTTCAGACATCCAGAATTCCAAAAAATAGTATGAAAATTGTTCAATTAGCTAAGAAAAAACTAATCAGCAATAAgc                                                                                                                                                   >  1:390423/1‑149 (MQ=255)
      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:440641/150‑1 (MQ=255)
                                                                          ttcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGt                                                                            <  1:454445/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:1124859/1‑150 (MQ=255)
                                                                                                                                                   aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg  >  1:740862/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGG  >  NC_000913/987271‑987567

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: