Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,076,604 G→A 5.7% L153L (CTG→CTA) ‡ yihY → putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,076,6040GA 5.7% 572.2 / 14.0 176L153L (CTG→CTA) ‡yihYputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (83/83);  new base A (4/6);  total (87/89)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 186 total aligned reads displayed.

GTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAT  >  NC_000913/4076455‑4076753
                                                                                                                                                     |                                                                                                                                                     
  tGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGc                                                                                                                                                     <  2:370171/150‑1 (MQ=255)
   ggggCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCa                                                                                                                                                    >  1:1408232/1‑150 (MQ=255)
     ggCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGTGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAgg                                                                                                                                                  >  2:1741/1‑150 (MQ=255)
     ggCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAgg                                                                                                                                                  >  2:1027561/1‑150 (MQ=255)
         tGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATATCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGcc                                                                                                                                              <  2:282815/150‑1 (MQ=255)
         tGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGcc                                                                                                                                              <  1:1378038/150‑1 (MQ=255)
           cGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGcccc                                                                                                                                            <  1:711988/150‑3 (MQ=255)
            gggCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGcccca                                                                                                                                           <  2:160281/150‑4 (MQ=255)
               cTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGccccagcg                                                                                                                                        <  1:1017122/150‑7 (MQ=255)
               cTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCCCTGCCAGCAGCGGccccagcg                                                                                                                                        <  1:1402345/150‑7 (MQ=255)
                tGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTgg                                                                                                                                       >  1:1277409/1‑150 (MQ=255)
                 gATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGc                                                                                                                                      >  1:1418341/1‑150 (MQ=255)
                 gATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTCGCGCAGTAAACGAGCGCGACCAAAAATTTCCTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGc                                                                                                                                      <  2:1156283/150‑1 (MQ=255)
                  aTCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCg                                                                                                                                     <  2:249251/150‑1 (MQ=255)
                     gTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTGGTTTGCCGTGTAGTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATc                                                                                                                                  <  2:1109307/150‑1 (MQ=255)
                      tCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCa                                                                                                                                 <  1:171385/150‑1 (MQ=255)
                       cACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAg                                                                                                                                >  2:965655/1‑150 (MQ=255)
                        aCCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGt                                                                                                                               <  2:914326/150‑1 (MQ=255)
                         ccGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGtt                                                                                                                              >  1:649931/1‑150 (MQ=255)
                            cgTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCt                                                                                                                           <  2:612473/150‑1 (MQ=255)
                             gTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTa                                                                                                                          <  2:309167/150‑1 (MQ=255)
                               tATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTAtt                                                                                                                        <  2:1276170/150‑1 (MQ=255)
                                  tgttgATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGc                                                                                                                     <  1:1214709/150‑1 (MQ=255)
                                  tgttgATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGc                                                                                                                     <  1:491897/150‑1 (MQ=255)
                                     tgATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGctct                                                                                                                  >  2:1425861/1‑150 (MQ=255)
                                      gATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCtctc                                                                                                                 >  1:69185/1‑150 (MQ=255)
                                      gATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCtctc                                                                                                                 >  1:457163/1‑150 (MQ=255)
                                        tGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTccc                                                                                                               <  2:748450/150‑1 (MQ=255)
                                         gTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCt                                                                                                              <  2:1331370/150‑1 (MQ=255)
                                           aCTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTgc                                                                                                            <  1:1431370/150‑1 (MQ=255)
                                           aCTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTgc                                                                                                            >  1:209692/1‑150 (MQ=255)
                                           aCTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTgc                                                                                                            <  2:33603/150‑1 (MQ=255)
                                            cTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTgcg                                                                                                           >  2:132629/1‑150 (MQ=255)
                                            cTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTgcg                                                                                                           <  2:1074841/150‑1 (MQ=255)
                                             tCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTga                                                                                                                     >  2:429111/1‑138 (MQ=255)
                                             tCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGc                                                                                                                     <  1:429111/139‑1 (MQ=255)
                                              ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCt                                                                                                         <  2:1092329/150‑1 (MQ=255)
                                               cATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTg                                                                                                        >  2:860185/1‑150 (MQ=255)
                                                aTCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTgg                                                                                                       <  1:1070783/150‑1 (MQ=255)
                                                 tCGATTGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTggg                                                                                                      <  1:1193671/150‑1 (MQ=255)
                                                   gATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcg                                                                                                    <  1:1353852/150‑1 (MQ=255)
                                                   gATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcg                                                                                                    <  2:836578/150‑1 (MQ=255)
                                                     tAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcgag                                                                                                  <  1:739959/150‑1 (MQ=255)
                                                        cgcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgcgagcga                                                                                               <  1:1329146/150‑1 (MQ=255)
                                                          cgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAtc                                                                                             <  2:781000/150‑1 (MQ=255)
                                                            ttGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAtctc                                                                                           >  2:128184/1‑150 (MQ=255)
                                                              gAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCaa                                                                                         <  1:629441/150‑1 (MQ=255)
                                                                aTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATa                                                                                       >  2:881506/1‑150 (MQ=255)
                                                                 tACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATAc                                                                                      <  2:573647/150‑1 (MQ=255)
                                                                     aTCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTc                                                                                  <  2:69185/150‑1 (MQ=255)
                                                                        tGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATc                                                                               >  1:399001/1‑150 (MQ=255)
                                                                         ggCGCAGTAAACGAGCGCGACCAAAAATTTACTCGATCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGTGCTGGGCGAGCGATCTCAATACTGTCATCg                                                                              >  2:840477/1‑150 (MQ=255)
                                                                            gcAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACa                                                                           <  1:12460/150‑1 (MQ=255)
                                                                               gTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACg                                                                        >  2:508663/1‑150 (MQ=255)
                                                                                tAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGt                                                                       >  2:234059/1‑150 (MQ=255)
                                                                                tAAACGAGCCCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGATGGCGCCGATGATGGCAGTGGCCAGTCTGGAGATCAGTTCCTATTTGATCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGt                                                                       >  2:705309/1‑150 (MQ=255)
                                                                                  aaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgc                                                                     >  1:26064/1‑150 (MQ=255)
                                                                                  aaCGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgc                                                                     >  1:443847/1‑150 (MQ=255)
                                                                                    cGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctg                                                                   >  1:254446/1‑150 (MQ=255)
                                                                                     gAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctgc                                                                  <  2:1301407/150‑1 (MQ=255)
                                                                                         gcgACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAt                                                              >  1:616576/1‑150 (MQ=255)
                                                                                             ccAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAttttt                                                          >  2:1173773/1‑150 (MQ=255)
                                                                                             ccAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAttttt                                                          >  1:605433/1‑150 (MQ=255)
                                                                                              cAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTg                                                                                    <  1:1305255/123‑1 (MQ=255)
                                                                                              cAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTc                                                         >  2:1237262/1‑150 (MQ=255)
                                                                                              cAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTc                                                         >  2:265456/1‑150 (MQ=255)
                                                                                               aaaaaTTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTcc                                                        <  1:132629/150‑1 (MQ=255)
                                                                                                ataaTTTACTCGTTCGACGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCg                                                       <  1:891052/148‑1 (MQ=255)
                                                                                                aaaaTTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCg                                                       >  1:766097/1‑150 (MQ=255)
                                                                                                  aaTTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGtt                                                     >  2:1075448/1‑150 (MQ=255)
                                                                                                  aaTTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGtt                                                     <  2:603711/150‑1 (MQ=255)
                                                                                                      tACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCATTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtg                                                 <  1:1161717/150‑1 (MQ=255)
                                                                                                       aCTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgt                                                <  1:1027561/150‑1 (MQ=255)
                                                                                                        cTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctgc                                                                  <  1:496190/131‑1 (MQ=255)
                                                                                                        cTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgctgc                                                                  >  2:496190/1‑131 (MQ=255)
                                                                                                        cTCGTTCGCCGTGTACTGGATGATTTTAACGCTGCGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgtt                                               >  2:982128/1‑150 (MQ=255)
                                                                                                          cGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgttgt                                             <  2:399001/150‑1 (MQ=255)
                                                                                                           gTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTc                                            <  2:254446/150‑1 (MQ=255)
                                                                                                           gTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTc                                            >  1:854487/1‑150 (MQ=255)
                                                                                                              cGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTg                                         <  2:671840/150‑1 (MQ=255)
                                                                                                                      aCTGGATGATTTTAAGGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCt                                 <  2:649931/150‑1 (MQ=255)
                                                                                                                          gatgatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTg                             >  1:800854/1‑150 (MQ=255)
                                                                                                                            tgatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGt                           <  2:531417/150‑1 (MQ=255)
                                                                                                                             gatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGtt                          <  2:800854/150‑1 (MQ=255)
                                                                                                                               ttttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTTCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGc                        >  1:999844/1‑150 (MQ=255)
                                                                                                                                tttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCt                       >  2:626536/1‑150 (MQ=255)
                                                                                                                                 ttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTg                      <  2:1050650/150‑1 (MQ=255)
                                                                                                                                 ttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTg                      <  2:844494/150‑1 (MQ=255)
                                                                                                                                  tAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGt                     >  1:510247/1‑150 (MQ=255)
                                                                                                                                    aCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTAc                   >  1:965780/1‑150 (MQ=255)
                                                                                                                                       cTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGc                <  1:367479/150‑1 (MQ=255)
                                                                                                                                         ggggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCAt              >  2:1326974/1‑150 (MQ=255)
                                                                                                                                           ggCCCCTGCCAGCAGCGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTg            >  1:418105/1‑150 (MQ=255)
                                                                                                                                           ggCCCCTGCCAGCAGCGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTg            >  1:293039/1‑150 (MQ=255)
                                                                                                                                             ccGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGtt          >  1:275423/1‑150 (MQ=255)
                                                                                                                                              cGCTTCTGGCAGGGGCCAGTGTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGATCGATCTGAATACTGTCATCGACAACGTGATGCGTATTTTTCCGTTGGTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGGATTGTTc         <  2:1413865/150‑1 (MQ=255)
                                                                                                                                              cGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTc         >  1:793888/1‑150 (MQ=255)
                                                                                                                                                 tgctgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTa      <  1:1291575/150‑1 (MQ=255)
                                                                                                                                                  gctgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAc     >  2:439636/1‑150 (MQ=255)
                                                                                                                                                     ggCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAt  >  1:1142763/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGGGGCGTGCGGGCTGATCGTCACCGCGTTATTGTTGATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACCAT  >  NC_000913/4076455‑4076753

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: