Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,440,356 A→T 5.0% L287Q (CTG→CAG)  ytfL ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,440,3560AT 5.0% 297.9 / 3.0 94L287Q (CTG→CAG) ytfLputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (37/52);  new base T (2/3);  total (39/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGATCGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCGGATGTTCCGCCA  >  NC_000913/4440207‑4440503
                                                                                                                                                     |                                                                                                                                                   
ccAGATCGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCTGCACAATCa                                                                                                                                                     <  2:751046/150‑1 (MQ=255)
ccAGATCGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCa                                                                                                                                                     <  1:735422/150‑1 (MQ=255)
  aGATCGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGc                                                                                                                                                   <  1:657539/150‑1 (MQ=255)
   gATCGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCACCGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCg                                                                                                                                                  >  2:116621/1‑150 (MQ=255)
      cGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGTGTTAACGTATCCGGCACAATCAGCgtgt                                                                                                                                               <  2:922928/150‑1 (MQ=255)
      cGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCgtgt                                                                                                                                               <  1:149147/150‑1 (MQ=255)
      cGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCgtgt                                                                                                                                               <  1:143352/150‑1 (MQ=255)
         ccATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGc                                                                                                                                            <  1:652519/150‑1 (MQ=255)
          cATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCg                                                                                                                                           >  1:642149/1‑150 (MQ=255)
           aTCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGa                                                                                                                                          <  1:529065/150‑1 (MQ=255)
                cGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTg                                                                                                                                     <  2:775913/150‑1 (MQ=255)
                  tGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTa                                                                                                                                   <  2:1184076/150‑1 (MQ=255)
                  tGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTa                                                                                                                                   >  1:1210388/1‑150 (MQ=255)
                  tGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTa                                                                                                                                   >  1:1204484/1‑150 (MQ=255)
                     tcatcaCCTCATTGAGGGTGATGATCACCCACACCAGCGCGTACTCGTTCAACATCACAGAGAAGTCTTCACCTGCGGTATTAAAAATTTCCAACGCCTCTGAAAGCGTTAACCTATCCGGCACAATAAGCGTGATGAGAATAAGTAcgc                                                                                                                                >  2:715035/1‑150 (MQ=255)
                      catcaCGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTAcgcc                                                                                                                               <  2:193405/150‑1 (MQ=255)
                      catcaCGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTAcgcc                                                                                                                               >  1:734038/1‑150 (MQ=255)
                        tcaCGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTAcgccgc                                                                                                                             <  2:1204619/150‑1 (MQ=255)
                        tcaCGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTAcgccgc                                                                                                                             >  1:520804/1‑150 (MQ=255)
                          aCGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCCCCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTg                                                                                                                           <  1:1370820/150‑1 (MQ=255)
                                ttGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGt                                                                                                                     <  1:1125085/150‑1 (MQ=255)
                                 tGAGGGTGATGATCCCCACCACCAGCGCGGACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGGGAATTTGTACGCCGCTGTTCAGTg                                                                                                                    <  2:127678/150‑1 (MQ=255)
                                   aGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGcc                                                                                                                  >  1:1317696/1‑150 (MQ=255)
                                     ggTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAATGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAg                                                                                                                <  2:938038/150‑1 (MQ=255)
                                     ggTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAg                                                                                                                <  2:504827/150‑1 (MQ=255)
                                      gTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCTGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAgg                                                                                                               >  1:197578/1‑150 (MQ=255)
                                      gTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCTGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAgg                                                                                                               >  1:191917/1‑150 (MQ=255)
                                       tgatgaTCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGc                                                                                                              <  2:1166513/150‑1 (MQ=255)
                                           gatCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTg                                                                                                          >  2:77061/1‑150 (MQ=255)
                                            atCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTTCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGACGCTGTTCAGTGCCAGGCTTTgg                                                                                                         >  2:649792/1‑150 (MQ=255)
                                            atCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTgg                                                                                                         >  1:1098891/1‑150 (MQ=255)
                                             tCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACTATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGt                                                                                                        <  2:61505/150‑1 (MQ=255)
                                             tCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGt                                                                                                        <  2:328853/150‑1 (MQ=255)
                                             tCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACCTATCCGGAACAATCAGCGTGTTGCGAATTTGTACGCCGCTCTTCAGTGCCAGGCTTTGGt                                                                                                        >  2:326403/1‑150 (MQ=255)
                                              ccccACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGtt                                                                                                       <  2:383146/150‑1 (MQ=255)
                                               cccACCACCAGCGCGTACTAGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTGGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTa                                                                                                      >  1:1111986/1‑150 (MQ=255)
                                                  accaccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGcc                                                                                                   <  1:585007/150‑1 (MQ=255)
                                                  accaccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGcc                                                                                                   <  1:179181/150‑1 (MQ=255)
                                                    caccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTa                                                                                                      >  1:610565/1‑145 (MQ=255)
                                                    caccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTa                                                                                                      <  2:610565/145‑1 (MQ=255)
                                                    caccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAg                                                                                                 >  1:298548/1‑150 (MQ=255)
                                                    caccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAg                                                                                                 <  2:618376/150‑1 (MQ=255)
                                                    caccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAg                                                                                                 <  2:881405/150‑1 (MQ=255)
                                                    caccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAg                                                                                                 >  1:681606/1‑150 (MQ=255)
                                                     accaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGGTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGc                                                                                                <  1:1177151/150‑1 (MQ=255)
                                                      ccaGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCa                                                                                               <  2:1505947/150‑1 (MQ=255)
                                                          cgcgTACTCGTTCATGATCACCGCGTAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGGGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCAcgcg                                                                                           <  2:1438476/150‑1 (MQ=255)
                                                          cgcgTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCAcgcg                                                                                           >  1:493674/1‑150 (MQ=255)
                                                               aCTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTca                                                                                      <  2:296747/150‑1 (MQ=255)
                                                                cTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTcag                                                                                     >  2:1084054/1‑150 (MQ=255)
                                                                 tCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTcagc                                                                                    <  2:937604/150‑1 (MQ=255)
                                                                  cGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTcagca                                                                                   <  1:183441/150‑1 (MQ=255)
                                                                      cATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTc                                                                               <  2:642149/150‑1 (MQ=255)
                                                                         gATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCttt                                                                            >  2:816870/1‑150 (MQ=255)
                                                                          aTCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTa                                                                           <  1:423108/150‑1 (MQ=255)
                                                                           tCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAg                                                                          >  2:41098/1‑150 (MQ=255)
                                                                            cACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGa                                                                         <  2:1221027/150‑1 (MQ=255)
                                                                               cgcgAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATc                                                                      <  2:520804/150‑1 (MQ=255)
                                                                               cgcgAAGTCTTCACCTGCGGTTTTAAAACTATCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATc                                                                      <  1:858091/150‑1 (MQ=255)
                                                                                gcgAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCGTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCGGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCg                                                                     <  2:411958/150‑1 (MQ=255)
                                                                                   aaGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACa                                                                  <  1:1116303/150‑1 (MQ=255)
                                                                                   aaGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACa                                                                  >  2:260603/1‑150 (MQ=255)
                                                                                       cTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAAc                                                              >  2:362827/1‑150 (MQ=255)
                                                                                            ccTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCgatg                                                         >  1:1024477/1‑150 (MQ=255)
                                                                                            ccTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCgatg                                                         <  2:790404/150‑1 (MQ=255)
                                                                                             cTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCgatga                                                        >  2:447552/1‑150 (MQ=255)
                                                                                              tGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCgatgat                                                       <  1:1591083/150‑1 (MQ=255)
                                                                                                   ttttAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGt                                                  >  1:492611/1‑150 (MQ=255)
                                                                                                    tttAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTc                                                 <  2:1617387/150‑1 (MQ=255)
                                                                                                    tttAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTc                                                 <  2:734038/150‑1 (MQ=255)
                                                                                                     ttAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCa                                                >  1:1508737/1‑150 (MQ=255)
                                                                                                           cTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCt                                          <  2:1317696/150‑1 (MQ=255)
                                                                                                                cAACGCCTCTTAAAGGGTTAACGTATCCTGCACAATCTGCTTGTTGCTAATTTGTACGCCGCTTTTCAGTGCCAGGCTTTGGTTAGCCAGCAGGCGGTTCAGCAGGTCTTTTGAATCGTCATAACCGATGATGTGGTCAATATCTTCAtt                                     <  1:715035/150‑1 (MQ=255)
                                                                                                                  aCGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTAc                                   <  1:77061/150‑1 (MQ=255)
                                                                                                                  aCGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTAc                                   <  1:1440673/150‑1 (MQ=255)
                                                                                                                       tctGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACg                              >  1:1455892/1‑150 (MQ=255)
                                                                                                                        ctGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGa                             >  2:1257515/1‑150 (MQ=255)
                                                                                                                        ctGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGa                             <  2:1098891/150‑1 (MQ=255)
                                                                                                                         ttaaaGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTAAGCAGGTCTTTAGAATCGACATAACCGATGAGGTGGTCAATATCTTCATGACAGACGAg                            >  2:976396/3‑150 (MQ=255)
                                                                                                                          gAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGGTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAgg                           >  2:221600/1‑150 (MQ=255)
                                                                                                                          gAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAgg                           <  1:462197/150‑1 (MQ=255)
                                                                                                                          gAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAgg                           >  2:1139973/1‑150 (MQ=255)
                                                                                                                                gTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATaa                                                               <  1:1052179/108‑1 (MQ=255)
                                                                                                                                gTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATaa                                                               >  2:1052179/1‑108 (MQ=255)
                                                                                                                                   aaCGTATCCGGCACAATCTGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTagag                  <  2:197578/150‑1 (MQ=255)
                                                                                                                                   aaCGTATCCGGCACAATCTGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTagag                  <  2:191917/150‑1 (MQ=255)
                                                                                                                                   aaCGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTagag                  >  1:1556654/1‑150 (MQ=255)
                                                                                                                                    aCGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGt                 >  1:1495348/1‑150 (MQ=255)
                                                                                                                                       tATCCGGCACAATCAGCGTGTTGCGAATTTGTACTCCGCTGTTCAGTGCCAGGGTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACGGATTATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCg              <  2:1266263/150‑1 (MQ=255)
                                                                                                                                       tATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCg              >  1:842682/1‑150 (MQ=255)
                                                                                                                                        aTCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGGTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCgg             <  2:1077910/150‑1 (MQ=255)
                                                                                                                                            ggCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCgg             <  1:936516/146‑1 (MQ=255)
                                                                                                                                            ggCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCgg             >  2:936516/1‑146 (MQ=255)
                                                                                                                                              cacaATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCGGATGTTc       >  2:919748/1‑150 (MQ=255)
                                                                                                                                                caATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCGGATGTTCCg     <  1:362827/150‑1 (MQ=255)
                                                                                                                                                   tCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCGGATGTTCCGCCa  <  1:447552/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCAGATCGCCCATCAGCGTGGTCATCACGTCATTGAGGGTGATGATCCCCACCACCAGCGCGTACTCGTTCATGATCACCGCGAAGTCTTCACCTGCGGTTTTAAAACTTTCCAACGCCTCTGAAAGGGTTAACGTATCCGGCACAATCAGCGTGTTGCGAATTTGTACGCCGCTGTTCAGTGCCAGGCTTTGGTTAGCCAGCACGCGGTTCAGCAGGTCTTTAGAATCGACATAACCGATGATGTGGTCAATATCTTCATTACAGACGAGGAACTTAGAGTGCGGATGTTCCGCCA  >  NC_000913/4440207‑4440503

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