Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,881,596 A→C 6.1% S183A (TCC→GCC) ‡ casC ← CRISP RNA (crRNA) containing Cascade antiviral complex protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,881,5960AC 6.1% 330.2 / 10.5 115S183A (TCC→GCC) ‡casCCRISP RNA (crRNA) containing Cascade antiviral complex protein
Reads supporting (aligned to +/- strand):  ref base A (49/59);  new base C (2/5);  total (51/64)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 115 total aligned reads displayed.

AGGTTAATGTTGGCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCTTTTTATCATCCAGA  >  NC_000913/2881448‑2881744
                                                                                                                                                    |                                                                                                                                                    
aGGTTAATGTTGGCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGAc                                                                                                                                                     >  1:660264/1‑150 (MQ=255)
 ggTTAATGTTGGCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACa                                                                                                                                                    >  2:1100532/1‑150 (MQ=255)
    tAATGTTGGCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTg                                                                                                                                                 >  2:64906/1‑150 (MQ=255)
       tGTTGGCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCAc                                                                                                                                              <  1:25916/150‑1 (MQ=255)
        gTTGGCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCAcc                                                                                                                                             <  1:534404/150‑1 (MQ=255)
            gCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGTGCAATGTCCATTGcgcatgcg                                                                                                                                         <  1:634495/150‑8 (MQ=255)
            gCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGTGCAATGTCCATTGcgcatgcg                                                                                                                                         <  1:1235116/150‑8 (MQ=255)
            gCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGTCCATTGcgcatgcg                                                                                                                                         <  2:156675/150‑8 (MQ=255)
                 aCGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGTCCATTGcgcatg                                                                                                                                           >  2:1203071/1‑138 (MQ=255)
                   gATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTc                                                                                                                                  <  2:1573745/150‑1 (MQ=255)
                    atatAAAACACCCGATGATAATTCCTGAGTTCCCAGATGTGCAGTACCTTGTTCCTGTAAATCATCTACAGCGGTGACCCAGTCAATATCAGAATCAACCTGATGAGTAGTTATCGCATGCGCTATGGACATTGCACAATCAACTTTTcc                                                                                                                                 <  2:416129/150‑1 (MQ=255)
                     tAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTccc                                                                                                                                <  2:1601233/150‑1 (MQ=255)
                     tAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTccc                                                                                                                                <  1:1468766/150‑1 (MQ=255)
                      aaaaaaCACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCa                                                                                                                               <  2:980686/150‑1 (MQ=255)
                        aaaaCACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAAc                                                                                                                             >  1:1257546/1‑150 (MQ=255)
                           acacCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCa                                                                                                                          >  2:117030/1‑150 (MQ=255)
                            cacCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCGCCATCAACTTTTCCCAACTCAg                                                                                                                         <  2:1542998/150‑1 (MQ=255)
                             acCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCACCTCAGt                                                                                                                        >  2:156477/1‑150 (MQ=255)
                              cccGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGtc                                                                                                                       >  2:128479/1‑150 (MQ=255)
                              cccGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGtc                                                                                                                       <  1:1633424/150‑1 (MQ=255)
                               ccGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCATAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGtca                                                                                                                      <  2:938488/150‑1 (MQ=255)
                               ccGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTAAAGCGGTGAACCAGTCAATATCAGACTCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGtca                                                                                                                      <  2:932357/150‑1 (MQ=255)
                               ccGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTACCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCATAATCAACCTGATGAGTTGTGATCGCATGCGCAATGGACATTGCACCAGCAACTTTTCCCAACTCAGtca                                                                                                                      >  2:318095/1‑150 (MQ=255)
                                cGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGtcat                                                                                                                     >  1:1186862/1‑150 (MQ=255)
                                cGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATAACATGCGCAATGCACATATCACCATCAACATTTCCCAACTCAGtcat                                                                                                                     >  2:1044377/1‑150 (MQ=255)
                                 gATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGtcatc                                                                                                                    >  2:322500/1‑150 (MQ=255)
                                  aTGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGtcatca                                                                                                                   <  1:343816/150‑1 (MQ=255)
                                   tGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTcatcat                                                                                                                  >  1:159169/1‑150 (MQ=255)
                                      aaaTTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGcc                                                                                                               <  2:232126/150‑1 (MQ=255)
                                        aTTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGc                                                                                                             >  2:1425517/1‑150 (MQ=255)
                                        aTTCCTGAGTTCCAGATGTGCAGAACCTTGTTACTGTAAATAATCTACCGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCt                                                                                                            >  2:906573/1‑150 (MQ=255)
                                            cTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGt                                                                                                         <  2:1645711/150‑1 (MQ=255)
                                             tGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGtt                                                                                                        >  1:676078/1‑150 (MQ=255)
                                                gTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCTGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGcc                                                                                                     <  1:1592894/150‑1 (MQ=255)
                                                gTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGcc                                                                                                     >  2:72603/1‑150 (MQ=255)
                                                 ttCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCa                                                                                                    <  2:660264/150‑1 (MQ=255)
                                                     cAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCAAAGTCAGTCATCATGCCGCTGGTTGACCttct                                                                                                >  2:173676/1‑150 (MQ=255)
                                                           tgCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACt                                                                                          <  2:1186862/150‑1 (MQ=255)
                                                           tgCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACt                                                                                          >  1:1461982/1‑150 (MQ=255)
                                                            gCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCATCTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTa                                                                                         >  2:1332613/1‑150 (MQ=255)
                                                              aGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAg                                                                                       >  2:1357489/1‑150 (MQ=255)
                                                              aGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAg                                                                                       <  2:1110689/150‑1 (MQ=255)
                                                               gAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAgc                                                                                      >  1:697879/1‑150 (MQ=255)
                                                                 aCCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAgcgc                                                                                    >  2:1402191/1‑150 (MQ=255)
                                                                  ccTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCa                                                                                   >  1:530197/1‑150 (MQ=255)
                                                                  ccTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCa                                                                                   >  1:1590202/1‑150 (MQ=255)
                                                                   cTTGTTCCTGTTAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCaa                                                                                  <  1:212419/150‑1 (MQ=255)
                                                                    ttGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAat                                                                                 <  1:737331/150‑1 (MQ=255)
                                                                     tGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCaaaa                                                                                >  1:1214859/1‑148 (MQ=255)
                                                                      tttCCTGTAAATCATCTACATCGGTGATCCAGTGAATATCAGAATCAACCTGATGAGTTGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGAAatat                                                                               <  1:1247765/149‑1 (MQ=255)
                                                                       ttCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATc                                                                              <  2:255066/150‑1 (MQ=255)
                                                                       ttCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATc                                                                              <  2:87167/150‑1 (MQ=255)
                                                                        tCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACATTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCa                                                                             <  1:150155/150‑1 (MQ=255)
                                                                          cTGTAAATCATCTACAGCGGTGCACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCAGCATGCCGCTGGTTGCCATTCTTCCACGTCTCGCAATATCAac                                                                           >  2:105000/1‑150 (MQ=255)
                                                                                atcatcTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCTATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCctg                                                                     <  1:1433913/150‑1 (MQ=255)
                                                                                atcatcTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGc                                                                                                      <  2:1034975/117‑1 (MQ=255)
                                                                                atcatcTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGc                                                                                                      >  1:1034975/1‑117 (MQ=255)
                                                                                 tcatcTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCctgc                                                                    <  1:1108030/150‑1 (MQ=255)
                                                                                     cTACAGCGGTGAACCAGTCAATATCAGAATCATCCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATTCCGCTGGTTGCCATTCTGCCACTAAGCGCAATATCAACACCCTGCTGTa                                                                <  1:1165474/150‑1 (MQ=255)
                                                                                     cTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTa                                                                >  1:776347/1‑150 (MQ=255)
                                                                                     atACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTa                                                                <  1:1663397/149‑1 (MQ=255)
                                                                                         aGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACTCCCTGCTGTAAAtt                                                            <  2:1539890/150‑1 (MQ=255)
                                                                                          gCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTc                                                           <  2:1214859/150‑1 (MQ=255)
                                                                                           cGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCAc                                                                                           >  2:1393657/1‑117 (MQ=255)
                                                                                           cGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCAc                                                                                           <  1:1393657/117‑1 (MQ=255)
                                                                                               gAAACAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTTCACAATCAACGTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCCGTACGCGCAATCTCAAAAAGCTGCTGGAAATTCACACg                                                      >  2:486011/1‑150 (MQ=255)
                                                                                                 aCCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACCCgga                                                    >  1:270985/1‑148 (MQ=255)
                                                                                                 aCCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAAAGGACATTGCACCAACAACTTTACCCAACTCAGTCATCATGCCGCTGCTTGCCATTCTTCCACTAAGCACAATATCAACACCCTGCAGTAAATTGACACGTa                                                    >  1:176640/1‑150 (MQ=255)
                                                                                                     ttCAATATCAGAATGAACCTGATGAGTAGTGTTCGCATGCGCAATGGAGATTGCACCATCATCTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCTCAATATCAACACCCTGCTGTAAATTCACACGTATggc                                                <  2:736401/149‑1 (MQ=255)
                                                                                                     gTCAATATCAGAAGCACCCTGATGAGTAGTGATCGCATGCGCAATGGACAGTGCACCTTCAACTTTTCACAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCGACAAAGCGCAATATCATCCCCGTGCAGTAAAGTCACCCGTAGggc                                                >  2:253510/1‑150 (MQ=255)
                                                                                                       cAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATggcgg                                              <  1:906573/150‑1 (MQ=255)
                                                                                                         atatCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCGATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCa                                            <  1:338005/150‑1 (MQ=255)
                                                                                                          tatCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCaa                                           >  2:1249599/1‑150 (MQ=255)
                                                                                                           atCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCCCACGTATTGCGGCAat                                          >  2:380069/1‑150 (MQ=255)
                                                                                                               gAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCt                                      <  1:1666036/150‑1 (MQ=255)
                                                                                                                 aTCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTc                                    <  2:1368870/150‑1 (MQ=255)
                                                                                                                 aTCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTc                                    <  1:1321992/150‑1 (MQ=255)
                                                                                                                  tCAACCTTATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTcc                                   <  2:590956/150‑1 (MQ=255)
                                                                                                                  tCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTcc                                   <  1:1702637/150‑1 (MQ=255)
                                                                                                                      ccTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTaa                               <  1:1402191/150‑1 (MQ=255)
                                                                                                                        tgatgaGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGa                             <  1:1327499/150‑1 (MQ=255)
                                                                                                                        tgatgaGTAGTGATAGCATGCGCAATTGACATTGGACCATCAACTTTTCCCAAATCAGTCATCATTCCTCTGGTTGCCATTCATCCACTAAGCGCAATATCAACACCCTTCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGa                             <  2:176640/150‑1 (MQ=255)
                                                                                                                         gatgaGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGaa                            >  2:1317336/1‑150 (MQ=255)
                                                                                                                            tagttgtGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTTTGGCGGCAATATCTTCCTTTAGAACtt                         <  1:1044377/149‑1 (MQ=255)
                                                                                                                                 gtGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGagc                    <  1:655904/150‑1 (MQ=255)
                                                                                                                                 gtGATCGCATGCGCAATGGACATTGAAACATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCATTATCAACACCCTGCTGTAAATTCACACGTATGGCGTCAATATCTTCCTTAAGAACTTTGagc                    <  1:268198/150‑1 (MQ=255)
                                                                                                                                   gATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGagcag                  >  1:1440030/1‑150 (MQ=255)
                                                                                                                                    aTCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGagcagc                 <  2:1461982/150‑1 (MQ=255)
                                                                                                                                    aTCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGagcagc                 <  1:117030/150‑1 (MQ=255)
                                                                                                                                      cGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCtt               >  2:645773/1‑150 (MQ=255)
                                                                                                                                      cGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCtt               <  1:1249599/150‑1 (MQ=255)
                                                                                                                                        cATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCtttt             >  2:844554/1‑150 (MQ=255)
                                                                                                                                        cATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCtttt             <  1:1100532/150‑1 (MQ=255)
                                                                                                                                         tggtgcAATGTCCATTGCGCCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGag                     >  2:402195/5‑141 (MQ=255)
                                                                                                                                         tggtgcAATGTCCATTGCGCCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGag                     <  1:402195/137‑1 (MQ=255)
                                                                                                                                              cAATGGACATTGCACCATCACCTTTTCCCAACTCACTCATCATGCCGCAGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAACTTCACACGTATGGCGGCAATTTCTTCCGAAAGAACTTTGAGCAGCTTTTTatcca       >  2:674272/1‑148 (MQ=255)
                                                                                                                                              cAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCTTTTTatcat       >  2:365462/1‑150 (MQ=255)
                                                                                                                                                  ggACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCTTTTTATCATCCAg   <  2:1590202/150‑1 (MQ=255)
                                                                                                                                                  ggACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCTTTTTATCATCCAg   >  1:42746/1‑150 (MQ=255)
                                                                                                                                                   gACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCTTTTTATCATCCAGa  >  1:1702938/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGTTAATGTTGGCATAACGATAAAAAACACCCGATGAAAATTCCTGAGTTCCCAGATGTGCAGAACCTTGTTCCTGTAAATCATCTACAGCGGTGAACCAGTCAATATCAGAATCAACCTGATGAGTAGTGATCGCATGCGCAATGGACATTGCACCATCAACTTTTCCCAACTCAGTCATCATGCCGCTGGTTGCCATTCTTCCACTAAGCGCAATATCAACACCCTGCTGTAAATTCACACGTATGGCGGCAATATCTTCCTTAAGAACTTTGAGCAGCTTTTTATCATCCAGA  >  NC_000913/2881448‑2881744

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