Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,770,319 C→G 16.5% P206R (CCG→CGG)  yfjR → CP4‑57 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,770,3190CG16.5% 300.8 / 51.2 115P206R (CCG→CGG) yfjRCP4‑57 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (53/43);  new base G (4/15);  total (57/58)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 119 total aligned reads displayed.

GGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTA  >  NC_000913/2770220‑2770418
                                                                                                   |                                                                                                   
  cATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCgg                                                                                                   <  1:2067018/100‑1 (MQ=255)
   atatCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGt                                                                                                  >  2:847610/1‑100 (MQ=255)
     atCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGATACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTCAAAAACCg                                                                                                    >  2:74342/1‑96 (MQ=255)
      tCGCCGTATTCACACTTGATGAAATCCATCTAATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCtt                                                                                               <  1:2548096/100‑1 (MQ=255)
        gCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttca                                                                                             >  2:786639/1‑100 (MQ=255)
         ccGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcat                                                                                            <  1:159980/100‑1 (MQ=255)
           gTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCGGGTCttcaacc                                                                                          <  1:2034440/100‑4 (MQ=255)
              ttCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAgg                                                                                       >  1:1256976/1‑100 (MQ=255)
              ttCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAgg                                                                                       <  2:2176120/100‑1 (MQ=38)
              ttCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATCACAGTCTGTGTCGCCTGGTTGAAGACCCTGTCttcattcatg                                                                                       <  2:1755171/100‑3 (MQ=255)
               tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCGGGTCTTCAACCAGGc                                                                                      <  1:748616/100‑1 (MQ=255)
               tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCttcaaccaggc                                                                                      <  2:1353787/100‑8 (MQ=37)
               tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCttcaaccaggc                                                                                      <  2:520170/100‑8 (MQ=37)
               tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTAttcaaccaggc                                                                                      <  1:1025648/100‑8 (MQ=37)
               tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAAGGAAGACCGGGTCttcaaccaggc                                                                                      <  2:2194998/100‑8 (MQ=38)
               tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGACTGAATGAAGACCGGGTCttcaaccaggc                                                                                      <  1:1046141/100‑8 (MQ=37)
                cacaCTTGATGAAATCCATCTGATTCATCCGTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGcg                                                                                     <  2:1229403/100‑2 (MQ=37)
                cacaCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                     <  1:1673634/100‑1 (MQ=255)
                cacaCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                     >  1:1984477/1‑100 (MQ=255)
                cacaCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                     >  1:2151882/1‑100 (MQ=255)
                cacaCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCttcaaccaggcg                                                                                     <  2:2126175/100‑9 (MQ=37)
                 acacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCttcaaccaggcga                                                                                    <  1:1642234/100‑10 (MQ=37)
                  cacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCtt                                                                                   >  1:2352973/1‑100 (MQ=255)
                  cacTTGATGAAATCCATCTGATTAAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCt                                                                                    >  2:1540163/1‑99 (MQ=255)
                     ttGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAcc                                                                                <  2:1426944/100‑1 (MQ=255)
                     ttGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGcgacaca                                                                                <  1:905753/100‑7 (MQ=25)
                      tgatgaGATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGAATGAAGACCGGGTCTTCAACCAGGcgacacag                                                                               <  1:393/100‑8 (MQ=25)
                      tgatgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAc                                                                                 >  1:1230182/1‑98 (MQ=255)
                      tgatgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTCTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAccc                                                                               >  1:1832491/1‑100 (MQ=255)
                      tgatgaAATCCATCTGATTCAGCCTTTGCAGGAGACGTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAccc                                                                               <  1:2203655/100‑1 (MQ=255)
                       gatgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAcccc                                                                              <  2:1996823/100‑1 (MQ=255)
                        atgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCa                                                                             <  2:1822213/100‑1 (MQ=255)
                         tgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCAt                                                                            >  2:1816752/1‑100 (MQ=255)
                               ccATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGc                                                                      >  2:1892592/1‑100 (MQ=255)
                                cATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                     <  1:275530/100‑1 (MQ=255)
                                cATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                     >  1:1295007/1‑100 (MQ=255)
                                cATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                     <  1:1545907/100‑1 (MQ=255)
                                cATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                     <  1:1867898/100‑1 (MQ=255)
                                cATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                     <  2:1132849/100‑1 (MQ=255)
                                cATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGATTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                     >  1:1822159/1‑100 (MQ=255)
                                cATCTGATTCAGCCTTTGCAGGAGACTTTGCGCTGCAATGACAGTCTGTGTCGCCAGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTAGc                                                                      >  2:232594/1‑99 (MQ=255)
                                 aTCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                    >  1:488160/1‑100 (MQ=255)
                                 aTCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCtt                                                                    >  2:569018/1‑100 (MQ=255)
                                  tCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCttt                                                                   >  2:90121/1‑100 (MQ=255)
                                   cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTTAAGACGCGGTCTTCATTCAGGGCTTCCCCCATTTTCGCTTTa                                                                  >  1:1158358/1‑100 (MQ=255)
                                   cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTa                                                                  >  2:1335809/1‑100 (MQ=255)
                                    tGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTAt                                                                 <  2:2451764/100‑1 (MQ=255)
                                      aTTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATcc                                                               <  2:265935/100‑1 (MQ=255)
                                         cAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGc                                                            >  1:1837146/1‑100 (MQ=255)
                                         cAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGc                                                            >  1:642614/1‑100 (MQ=255)
                                          aGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTACATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCa                                                           <  2:1820335/100‑1 (MQ=255)
                                           gCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCAt                                                          <  1:273720/100‑1 (MQ=255)
                                             cTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTc                                                        >  2:2137912/1‑100 (MQ=255)
                                             cTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTc                                                        >  1:2115051/1‑100 (MQ=255)
                                              tttGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGCCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCa                                                       <  2:1821384/100‑1 (MQ=255)
                                               ttGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCAc                                                      <  2:143283/100‑1 (MQ=255)
                                                tGCAGGAGACTTTTCTCCGTATTGACAGTCTGTGTCTCCTGGTTTAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCAc                                                      <  1:1498474/99‑1 (MQ=255)
                                                  cAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGc                                                   >  2:889697/1‑100 (MQ=255)
                                                  cAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGc                                                   >  2:889776/1‑100 (MQ=255)
                                                  cAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGc                                                   >  2:2476374/1‑100 (MQ=255)
                                                    ggAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCtt                                                 >  1:1636375/1‑100 (MQ=255)
                                                    gcagaCTTTTCGCCGCAATGACAGTCTGTCTCGCCTGGTTGAAGACCCGGTCTTCATTCAAGCCTTACCCCATTTTCGCTTTATCCAGCATTCCCTGCt                                                  >  1:991873/3‑99 (MQ=255)
                                                      agaCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTAc                                               <  2:1149587/100‑1 (MQ=255)
                                                         cTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGtt                                            >  2:2267440/1‑100 (MQ=255)
                                                         cTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGtt                                            <  2:833606/100‑1 (MQ=255)
                                                           tttCGCCGCAATGACAGTCTGTGTCGCCTTGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTg                                          <  1:1518787/100‑1 (MQ=255)
                                                           tttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTg                                          <  2:1636384/100‑1 (MQ=255)
                                                           tttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTg                                          >  2:1100238/1‑100 (MQ=255)
                                                              ngccgcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTc                                       >  1:259647/2‑100 (MQ=255)
                                                              cgccgcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTc                                       >  2:1883423/1‑100 (MQ=255)
                                                               gccgcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTcc                                      <  2:2422916/100‑1 (MQ=255)
                                                               gccgcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTcc                                      <  2:257590/100‑1 (MQ=255)
                                                                  gcAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGc                                   >  1:618960/1‑100 (MQ=255)
                                                                   cAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGcc                                  <  1:1518969/100‑1 (MQ=255)
                                                                     aTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTTTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCTCTGCTTACGTTTGTTCCGGCCg                                 <  1:735185/99‑1 (MQ=255)
                                                                      tGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGAc                               <  1:434867/100‑1 (MQ=255)
                                                                       gACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACa                              <  2:1498005/100‑1 (MQ=255)
                                                                        aCAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACan                             <  1:1872487/100‑2 (MQ=255)
                                                                        aCAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAg                             >  1:2448790/1‑100 (MQ=255)
                                                                         cAGTCTGTGTCGCCTGGTTGAAGCCCCCGTCTTCATTCAGGCCTTACCCCATTTTCGCTTGATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGc                            <  2:117522/100‑1 (MQ=255)
                                                                           gTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGccc                          >  2:1594606/1‑100 (MQ=255)
                                                                           gTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGccc                          <  2:2543305/100‑1 (MQ=255)
                                                                              tgtgTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCAcc                       <  2:583103/100‑1 (MQ=255)
                                                                              tgtgTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCAcc                       >  2:345065/1‑100 (MQ=255)
                                                                                  tCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGaa                   <  2:620688/100‑1 (MQ=255)
                                                                                   ggcctgaatGAAGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAAt                  >  2:1791657/9‑100 (MQ=255)
                                                                                    gcctgaatGAAGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATa                 >  2:376987/8‑100 (MQ=255)
                                                                                    gcctgaatGAAGACCGGGTCTTCAACCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATa                 >  2:1061599/8‑100 (MQ=37)
                                                                                     ccTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAg                >  2:1083898/1‑100 (MQ=255)
                                                                                      cTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgc               >  1:2444117/1‑100 (MQ=255)
                                                                                      cTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgc               >  1:1791409/1‑100 (MQ=255)
                                                                                      cTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgc               >  1:752990/1‑100 (MQ=255)
                                                                                        gaatGAAGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgcgc             >  2:1515095/4‑100 (MQ=255)
                                                                                           tGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAgg          >  2:1364147/1‑100 (MQ=255)
                                                                                             aaGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCg        <  2:493414/100‑1 (MQ=255)
                                                                                             aaGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCg        >  1:872217/1‑100 (MQ=255)
                                                                                              aGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTAACTGCTTACGTTTGTTTCGGCCGACAGCCCACCGGAATAGCGCAGGCgt       >  1:1003812/1‑100 (MQ=255)
                                                                                               gACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgtt      >  1:784085/1‑100 (MQ=255)
                                                                                                 cccGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACATCCCACCGGAATAGCGCAGGCgttgt    >  2:2027790/1‑100 (MQ=255)
                                                                                                 cccGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgt    <  2:1800449/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTA  >  NC_000913/2770220‑2770418

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