Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 573,991 A→G G12G (GGA→GGG quuD → DLP12 prophage; putative antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913573,9910AG98.5% 222.2 / ‑4.2 66G12G (GGA→GGGquuDDLP12 prophage; putative antitermination protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (45/20);  total (46/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CAAAACCCTGACAGGAGCAGGAGTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGAGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGCCCGTCTGA  >  NC_000913/573843‑574136
                                                                                                                                                    |                                                                                                                                                 
cAAAACCCTGACAGGAGCAGGAGTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGc                                                                                                                                                 >  2:663933/1‑151 (MQ=255)
cAAAACCCTGACAGGAGCAGGAGTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGc                                                                                                                                                 >  1:808967/1‑151 (MQ=255)
  aaaCCCTGACAGGAGCAGGAGTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGt                                                                                                                                               <  2:237936/151‑1 (MQ=255)
   aaCCCTGACAGGAGCAGGAGTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTg                                                                                                                                              >  1:682551/1‑151 (MQ=255)
     cccTGACAGGAGCAGGAGTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTggg                                                                                                                                            >  1:773528/1‑151 (MQ=255)
     cccTGACAGGAGCAGGAGTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTggg                                                                                                                                            >  2:623029/1‑151 (MQ=255)
            aggagcaggagTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCta                                                                                                                                      >  2:199248/1‑150 (MQ=255)
                 caggagTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATc                                                                                                                                >  1:224132/1‑151 (MQ=255)
                       tATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGa                                                                                                                            >  2:293509/1‑149 (MQ=255)
                        aTGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGa                                                                                                                            >  1:469367/1‑148 (MQ=255)
                         tGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGa                                                                                                                            >  1:496306/1‑147 (MQ=255)
                           aaGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGAtgt                                                                                                                      <  1:814748/151‑1 (MQ=255)
                               cAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGAcc                                                                                                                  >  2:401660/1‑151 (MQ=255)
                                agagATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCt                                                                                                                 >  2:344128/1‑151 (MQ=255)
                                      aaCAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGgtcgt                                                                                                           >  2:207943/1‑151 (MQ=255)
                                      aaCAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGgtcg                                                                                                            <  1:580719/150‑1 (MQ=255)
                                       aCAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGgtcgtc                                                                                                          >  2:19302/1‑151 (MQ=255)
                                       aCAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGgtcgtc                                                                                                          >  1:505248/1‑151 (MQ=255)
                                         aaTCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGAGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTcc                                                                                                         >  1:184737/1‑150 (MQ=255)
                                            cAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTTGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCAtt                                                                                                      <  2:148854/150‑1 (MQ=255)
                                             aaCTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTgc                                                                                                    >  1:90694/1‑151 (MQ=255)
                                             aaCTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTgc                                                                                                    >  2:151482/1‑151 (MQ=255)
                                              aCTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTgcc                                                                                                   <  2:447976/151‑1 (MQ=255)
                                              aCTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTgcc                                                                                                   >  2:549141/1‑151 (MQ=255)
                                               cTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTgccg                                                                                                  >  2:110677/1‑151 (MQ=255)
                                                 aaTACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTgccg                                                                                                  <  2:175751/149‑1 (MQ=255)
                                                    aCTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGt                                                                                             >  2:723635/1‑151 (MQ=255)
                                                        ttttGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGtttt                                                                                          >  2:448158/1‑150 (MQ=255)
                                                        ttttGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGtttt                                                                                          >  1:401213/1‑150 (MQ=255)
                                                         tttGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTa                                                                                         >  2:302947/1‑150 (MQ=255)
                                                          ttgttgATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTaa                                                                                        <  1:194365/150‑1 (MQ=255)
                                                           tgttgATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAg                                                                                       <  2:143861/150‑1 (MQ=255)
                                                            gttgATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAgg                                                                                      <  1:409733/150‑1 (MQ=255)
                                                            gttgATTTTTGCTTGTAATTGGCGTTCTGGGCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAggg                                                                                     >  2:77020/1‑151 (MQ=255)
                                                                aTTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGatt                                                                                  >  2:377829/1‑150 (MQ=255)
                                                                aTTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGatt                                                                                  >  2:459978/1‑150 (MQ=255)
                                                                   tttGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGattaatt                                                                              >  1:475988/1‑151 (MQ=255)
                                                                   tttGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGattaat                                                                               <  1:231597/150‑1 (MQ=255)
                                                                    ttGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGattaatt                                                                              >  1:829360/1‑150 (MQ=255)
                                                                        ttGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTc                                                                         >  2:774014/1‑151 (MQ=255)
                                                                          gTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCa                                                                        <  1:412544/150‑1 (MQ=255)
                                                                            aaTTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCaaa                                                                      >  2:317202/1‑150 (MQ=255)
                                                                            aaTTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCaaa                                                                      >  1:173823/1‑150 (MQ=255)
                                                                                       tGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCg                                                          >  2:35276/1‑151 (MQ=255)
                                                                                           cTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCg                                                      >  2:859778/1‑151 (MQ=255)
                                                                                           cTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCg                                                      >  2:506484/1‑151 (MQ=255)
                                                                                            tGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGc                                                     >  1:804740/1‑151 (MQ=255)
                                                                                            tGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGc                                                     >  1:34976/1‑151 (MQ=255)
                                                                                             gATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGc                                                     <  2:94940/150‑1 (MQ=255)
                                                                                                 tttGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGAGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAAtgt                                                >  1:60034/1‑151 (MQ=255)
                                                                                                    gtgGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAAtgttgt                                             <  1:360132/151‑1 (MQ=255)
                                                                                                            aGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATg                                      >  1:682025/1‑150 (MQ=255)
                                                                                                              ttGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGAcg                                   >  1:586277/1‑151 (MQ=255)
                                                                                                              ttGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGAc                                    <  2:586680/150‑1 (MQ=255)
                                                                                                                    cGTGATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCgatga                             <  2:808902/151‑1 (MQ=255)
                                                                                                                       gATATTCAGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATc                           <  1:14852/150‑1 (MQ=255)
                                                                                                                              aGATGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCgg                   <  1:754447/151‑1 (MQ=255)
                                                                                                                                aTGGTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCggg                  >  2:211945/1‑150 (MQ=255)
                                                                                                                                   gTTCTTGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGc               >  1:827302/1‑150 (MQ=255)
                                                                                                                                        tGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGcc         >  1:359444/1‑151 (MQ=255)
                                                                                                                                        tGAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGcc         >  1:227396/1‑151 (MQ=255)
                                                                                                                                         gAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGccc        >  1:222220/1‑151 (MQ=255)
                                                                                                                                         gAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGccc        <  2:847055/151‑1 (MQ=255)
                                                                                                                                         gAGCGTTGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGcc         <  2:496257/150‑1 (MQ=255)
                                                                                                                                              ttGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGCCCGTCTg   >  2:828042/1‑151 (MQ=255)
                                                                                                                                               tGGGGGGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGCCCGTCTg   <  2:717853/150‑1 (MQ=255)
                                                                                                                                                ggggggCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGCCCGTCTGa  <  1:623422/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAAAACCCTGACAGGAGCAGGAGTATGAAGTCAGAGATAACAATCAACTAATACTGTTTTGTTGATTTTTGCTTGTAATTGGCGTTCTGGTCTGATTTTTGTGGAGTAAGTTGATGCGTGATATTCAGATGGTTCTTGAGCGTTGGGGAGCGTGGGCGGCTAATAATCATGAAGATGTGACCTGGTCGTCCATTGCCGCCGGTTTTAAGGGATTAATTACTTCAAAAGTAAAATCTCGCCCGCAATGTTGTGACGATGACGCGATGATCATTTGCGGGTGCATGGCCCGTCTGA  >  NC_000913/573843‑574136

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