Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,775,089 T→A 5.2% intergenic (+68/‑456) yfjW → / → ypjI CP4‑57 prophage; putative inner membrane protein/pseudogene, CP4‑57 putative prophage remnant;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,775,0890TA 5.2% 280.9 / 3.8 76intergenic (+68/‑456)yfjW/ypjICP4‑57 prophage; putative inner membrane protein/pseudogene, CP4‑57 putative prophage remnant;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (40/32);  new base A (2/2);  total (42/34)
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

ATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGAA  >  NC_000913/2774970‑2775233
                                                                                                                       |                                                                                                                                                
attcattcattcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACt                                                                                                                     >  1:611598/1‑149 (MQ=255)
     ttcattcTTACTTTAAGAGTGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGCTTTAAAAACAGGc                                                                                                                                         <  4:777707/124‑1 (MQ=255)
     ttcattcTTACTTTAAGAGTGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGCTTTAAAAACAGGc                                                                                                                                         >  3:777707/1‑124 (MQ=255)
     ttcattcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGacac                                                                                                                      >  1:1112810/1‑143 (MQ=255)
     ttcattcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGacac                                                                                                                      <  2:1112810/143‑1 (MQ=255)
       cattcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCt                                                                                                              <  2:611598/149‑1 (MQ=255)
                tttAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGa                                                                                                                         <  4:1400842/129‑1 (MQ=255)
                tttAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGa                                                                                                                         >  3:1400842/1‑129 (MQ=255)
                     gAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTg                                                                                                                                <  2:380946/117‑1 (MQ=255)
                     gAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTg                                                                                                                                >  1:380946/1‑117 (MQ=255)
                                      aTTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCAt                                                                               >  2:945393/1‑149 (MQ=255)
                                                              tCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACt                                                       >  3:1271130/1‑149 (MQ=255)
                                                              tCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACt                                                       >  2:973808/1‑149 (MQ=255)
                                                              tCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTAAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACt                                                       >  4:1237355/1‑149 (MQ=255)
                                                                cAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGc                                                     <  1:973808/149‑1 (MQ=255)
                                                                cAATTGGAGCCTTTTTATGACTATTCATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGc                                                     <  3:1237355/149‑1 (MQ=255)
                                                                  aTTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTaaa                                                                                                                                               <  4:466531/57‑1 (MQ=255)
                                                                  aTTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTaaa                                                                                                                                               >  3:466531/1‑57 (MQ=255)
                                                                  aTTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCa                                                   >  3:1511981/1‑149 (MQ=255)
                                                                       aGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTaaa                                                                                                      >  1:214606/1‑93 (MQ=255)
                                                                       aGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTaaa                                                                                                      <  2:214606/93‑1 (MQ=255)
                                                                           tttttATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGttttt                                                                                           <  4:1044070/100‑1 (MQ=255)
                                                                           tttttATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGttttt                                                                                           >  3:1044070/1‑100 (MQ=255)
                                                                           tttttATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCtgt                                                                <  4:311154/127‑1 (MQ=255)
                                                                           tttttATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCtgt                                                                >  3:311154/1‑127 (MQ=255)
                                                                           tttttATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCt                                                                  <  3:569614/125‑1 (MQ=255)
                                                                           tttttATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCt                                                                  >  4:569614/1‑125 (MQ=255)
                                                                              ttATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGc                                                     >  4:162144/1‑135 (MQ=255)
                                                                              ttATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGc                                                     <  3:162144/135‑1 (MQ=255)
                                                                               tATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                      <  1:945393/149‑1 (MQ=255)
                                                                               tATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                      >  1:29901/1‑149 (MQ=255)
                                                                                  gACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTATAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                      <  4:1097591/146‑1 (MQ=255)
                                                                                  gACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                      >  3:1097591/1‑146 (MQ=255)
                                                                                    cTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaaaa                                 >  3:1775906/1‑149 (MQ=255)
                                                                                     tatttattAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAg                                <  2:29901/149‑1 (MQ=255)
                                                                                     tatttattAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAg                                >  2:980939/1‑149 (MQ=255)
                                                                                      atttattAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGc                               >  2:1879572/1‑149 (MQ=255)
                                                                                      atttattAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGc                               <  4:1271130/149‑1 (MQ=255)
                                                                                           ttaAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGt                                                                                                                                 >  4:1205206/1‑46 (MQ=255)
                                                                                           ttaAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGt                                                                                                                                 <  3:1205206/46‑1 (MQ=255)
                                                                                           ttaAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGaa                          >  1:1872815/1‑149 (MQ=255)
                                                                                           ttaAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGaa                          >  4:136233/1‑149 (MQ=255)
                                                                                            taAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATAc                                                                             >  3:1842362/1‑97 (MQ=255)
                                                                                            taAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATAc                                                                             <  4:1842362/97‑1 (MQ=255)
                                                                                            taAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAg                                                  >  2:1788048/1‑124 (MQ=255)
                                                                                            taAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAg                                                  <  1:1788048/124‑1 (MQ=255)
                                                                                              aaaTATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAg                                           >  2:288954/1‑129 (MQ=255)
                                                                                              aaaTATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAg                                           <  1:288954/129‑1 (MQ=255)
                                                                                              aaaTATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGt                                       <  2:1762816/133‑1 (MQ=255)
                                                                                              aaaTATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGt                                       >  1:1762816/1‑133 (MQ=255)
                                                                                              aaaTATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaaaa                                 >  2:1548660/1‑139 (MQ=255)
                                                                                              aaaTATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaaaa                                 <  1:1548660/139‑1 (MQ=255)
                                                                                                  atCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGTCAGAAATATGcc                   <  3:136233/149‑1 (MQ=255)
                                                                                                    ccTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTgttgtt                                                            >  2:1463571/1‑106 (MQ=255)
                                                                                                    ccTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTgttgtt                                                            <  1:1463571/106‑1 (MQ=255)
                                                                                                    ccTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTgttgtt                                                            <  3:634437/106‑1 (MQ=255)
                                                                                                    ccTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTgttgtt                                                            >  4:634437/1‑106 (MQ=255)
                                                                                                    ccTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGa                 <  4:1511981/149‑1 (MQ=255)
                                                                                                     cTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                      >  2:1040558/1‑127 (MQ=255)
                                                                                                     cTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                      <  1:1040558/127‑1 (MQ=255)
                                                                                                     cTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAatat                      >  4:504192/1‑143 (MQ=255)
                                                                                                     cTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAatat                      <  3:504192/143‑1 (MQ=255)
                                                                                                      tAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCgctg               >  4:1707443/1‑146 (MQ=255)
                                                                                                            aaGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCa                                                                                <  1:159737/78‑1 (MQ=255)
                                                                                                            aaGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCa                                                                                >  2:159737/1‑78 (MQ=255)
                                                                                                            aaGCCTGCTTTAAAAACAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCtat                                                                         <  2:1103095/85‑1 (MQ=39)
                                                                                                            aaGCCTGCTTTAAAAACAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCtat                                                                         >  1:1103095/1‑85 (MQ=39)
                                                                                                              gCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTg                                                                                                   <  1:58486/57‑1 (MQ=255)
                                                                                                              gCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTg                                                                                                   >  2:58486/1‑57 (MQ=255)
                                                                                                              gCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTc                                                                                        <  1:1059487/68‑1 (MQ=255)
                                                                                                              gCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTc                                                                                        >  2:1059487/1‑68 (MQ=255)
                                                                                                                cTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCa                                                   <  4:1025584/103‑1 (MQ=255)
                                                                                                                cTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCa                                                   >  3:1025584/1‑103 (MQ=255)
                                                                                                                cTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATg     >  3:1119666/1‑149 (MQ=255)
                                                                                                                 tGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATgg    >  4:1656391/1‑149 (MQ=255)
                                                                                                                 tGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATgg    >  4:493234/1‑149 (MQ=255)
                                                                                                                   tttttAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGaa  <  1:1879572/149‑1 (MQ=255)
                                                                                                                       |                                                                                                                                                
ATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTTAAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATACCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGAA  >  NC_000913/2774970‑2775233

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: